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Preface 


The  purpose  of  this  study  was  to  investigate  different 
techniques  for  selecting  the  optimum  downlink  signal  received 
by  spatially  diverse  earth  stations  operating  in  the  extreme¬ 
ly  high  frequency  (EHF )  band.  Space  diversity  configura¬ 
tions  can  be  used  to  overcome  the  attenuation  effects  of 
rainfall,  which  can  be  very  severe  for  frequencies  above 
10  GHz.  Thus,  special  consideration  must  be  given  to  the 
parameters  and  methods  used  to  determine  the  optimum  received 
signal  when  employing  a  space  diversity  reception  system. 

This  report  is  limited  in  scope  to  the  development  of  three 
linear  diversity  combining  techniques.  The  theoretical 
analysis  and  development  of  each  of  these  techniques  can  be 
fully  understood  by  the  reader  with  a  background  in  proba¬ 
bility  theory  and  digital  signaling  techniques. 

I  would  like  to  thank  Mr.  Richard  A.  Williams  of  the 
Defense  Communications  Engineering  Center  (DCEC)  for  pro¬ 
posing  this  project.  Gratitude  is  also  due  to  Captain  Jim 
Frazier,  also  from  DCEC,  for  his  assistance  in  arranging  our 
visits  to  DCEC.  Special  thanks  to  my  faculty  advisor.  Major 
Kenneth  G.  Castor  of  the  Air  Force  Institute  of  Technology, 
for  his  guidance  and  understanding  that  greatly  aided  me  in 
the  completion  of  this  study. 

Steven  L.  Skipper 
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Abstract 


This  study  investigated  some  of  the  diversity  combining 
techniques  that  can  be  used  with  a  space  diversity  downlink 
reception  configuration  to  overcome  the  attenuation  effects 
of  rainfall.  The  diversity  combining  techniques  selected 
for  development  and  analysis  were  selection,  equal  gain,  and 
maximal  ratio  combining.  Each  of  these  techniques  were  con¬ 
sidered  with  two  different  digital  signaling  schemes,  8-ary 
phase-shift  keying  (PSK)  and  differentially  encoded  quadra¬ 
ture  PSK  (DQPSK) ,  to  form  six  distinct  diversity  technique/ 
signaling  scheme  combinations. 

The  analysis  was  accomplished  by  deriving  the  probability 
of  bit  error  functions  for  8-ary  PSK  and  DQPSK  and  by  devel¬ 
oping  the  probability  of  density  functions  for  each  diversity 
technique  under  the  assumption  of  a  slowly  fading  Rayleigh 
channel.  By  taking  the  integral  from  zero  to  infinity  of  the 
product  of  the  two  functions  for  each  diversity  technique/ 
signaling  scheme  combination,  the  probability  of  bit  error 
was  calculated  for  each  specified  combination.  The  number  of 
channels  in  the  diversity  system  was  varied  to  illustrate  the 
trade-off  between  the  improvement  of  bit  error  performance 
and  cost.  The  results  show  the  theoretical  benefits  attrib¬ 
uted  to  the  diversity  combining  techniques  that  can  be 
employed  in  satellite  communications  systems.  Additionally, 
information  is  provided  to  estimate  the  rain  attenuation 


vii 


given  the  rain  rate  so  that  its  effects  on  bit  error  per 
formance  can  be  obtained  and  compared  to  the  case  of  no 
attenuation. 
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ANALYSIS  OF  SPACE  DIVERSITY  TECHNIQUES 


FOR  EHF  SATELLITE  DOWNLINKS 


I .  Introduction 


Background 

The  use  of  orbiting  satellites  has  become  an  integral 
part  of  our  world-wide  communications  and  military  command, 
control,  and  communications  systems.  A  satellite  communica¬ 
tions  system  can  be  categorized  as  one  of  three  basic  types 
shown  in  Figure  1.1.  System  I  shows  an  uplink  from  a  ground- 
based  earth  station  to  the  satellite,  and  a  downlink  from  the 
satellite  back  to  the  destination  earth  station.  The  func¬ 
tion  of  the  satellite  is  to  collect  the  electromagnetic 
fields  from  the  signal  transmitted  by  the  sending  earth  sta¬ 
tion  and  retransmit  the  received  signal  to  the  other  earth 
station . 

System  II  shows  a  satellite  crosslink  between  two  satel¬ 
lites  prior  to  the  downlink  transmission.  System  III  shows  a 
satellite  relay  system  involving  an  earth  station,  a  near- 
earth  user  (aircraft,  ship,  etc.)  and  a  satellite.  The  sys¬ 
tem  III  configuration  is  of  particular  importance  since  it 
has  many  military  applications.  Regardless  of  which  system 
configuration  is  employed,  there  are  several  atmospheric  ef¬ 
fects  that  can  degrade  the  transmitted  signal  strength  and 
reduce  the  capability  of  transmitting  information  via 


satellite. 


Figure  1.1  Types  of  Satellite  Systems 
I  Ground-ground 
II  Ground-crosslink -ground 
III  Ground-user  relay 

Source;  (11:5) 


The  most  serious  of  these  atmospheric  effects  to  a  sat¬ 
ellite  communications  link  is  rainfall,  which  affects  both 
uplink  and  downlink  transmissions  (Ili96).  Rain  effects 
become  more  severe  at  wavelengths  that  approach  the  water 
drop  size,  which  is  dependent  on  the  type  of  rainfall.  In 
heavy  rains  as  raindrop  size  increases,  severe  absorption 
occurs  in  the  extremely  high  frequency  (EHF)  band,  the  range 
where  many  existing  satellites  operate.  Therefore,  in  the 
EHF  band,  there  is  a  potential  for  significantly  more  signal 
attenuation  during  rain  as  compared  to  the  attenuation  caused 
during  clear  weather  conditions. 

Generally,  the  solution  to  maintaining  the  satellite 
communications  link  is  to  add  enough  extra  power  (power  mar¬ 
gin)  to  transmit  the  signal  over  the  maximum  additional 
attenuation  caused  by  rain.  However,  to  maintain  a  highly 
reliable  link  in  the  EHF  band,  the  cost  of  providing  suffi¬ 
cient  power  margin  with  current  technology  becomes  prohibi¬ 
tive  and  may  not  be  technically  feasible.  One  alternative 
for  improving  reliability  is  to  spatially  separate  two  or 
more  earth  stations  since  they  are  not  likely  to  experience 
simultaneous  fading  from  the  same  rainstorm  (See  Figure  1.2). 
For  the  uplink  transmission,  all  earth  station  -.ransmit  the 
same  signal  and  the  optimum  received  signal  (i.e.,  the  signal 
with  the  least  amount  of  attenuation)  is  chosen  at  the  satel¬ 
lite  and  retransmitted.  For  the  downlink  reception,  earth 
stations  receive  the  transmitted  signal  from  the  satellite 


simultaneously  and  the  optimum  signal  is  selected  and  then 
retransmitted  to  terrestrial  destinations. 

Problem 

When  using  a  space  diversity  configuration  to  overcome 
the  attenuation  effects  of  rainfall,  some  consideration  must 
be  given  to  the  parameters  and  the  methods  used  to  determine 
the  optimum  received  signal.  The  purpose  of  this  research 
is  to  investigate  techniques  for  selecting  the  optimum  signal 
from  signals  received  by  spatially  diverse  earth  stations 
operating  in  the  EHF  band  (i.e.,  the  case  of  Figure  1.2b). 

Scope 

The  scope  of  this  research  will  be  limited  to  an  inves¬ 
tigation  of  methods  to  determine  the  optimum  signal  received 
by  two  or  more  spatially  diverse  earth  stations  from  the  sat¬ 
ellite  downlink  transmissions.  Different  methods  of  selecting 
the  optimum  signal  will  be  investigated  and  developed.  After 
a  theoretical  analysis  of  each  selection  method,  a  computer 
program  will  be  developed  to  simulate  the  operational  environ¬ 
ment  and  verify  the  best  method  for  choosing  the  optimum  re¬ 
ceived  downlink  signal  in  a  space  diversity  configuration 
system.  Costs  associated  with  the  development  of  the  pro¬ 
posed  methods  for  signal  selection  will  not  be  examined. 

Assumptions 

The  following  assumptions  are  made  throughout  this 


research  project* 


1.  A  satellite  beacon  is  available. 


2.  All  uplink  transmissions  are  at  30  GHz  and  the  down¬ 
link  transmissions  are  at  20  GHz. 

3.  The  focus  of  this  report  will  be  on  satellite  sys¬ 
tems  using  diversity  techniques  to  improve  downlink  recep¬ 
tion  performance  as  depicted  in  Figure  1.2b. 

4.  Digital  signaling  schemes  used  will  be  8-ary  Phase- 
Shift  Keying  (PSK)  and  Differentially  encoded  Quadrature 

PSK  (DQPSK) . 

These  assumptions  are  established  to  model  a  generic 
satellite  system  for  military  applications  that  could  easily 
be  developed  and  deployed  using  current  satellite  system  tech¬ 
nology  and  practices . 

Summary  of  Current  Knowledge 

The  use  of  space  diversity  to  combat  atmospheric  attenua¬ 
tion  in  satellite  communications  links  is  a  well  known  concept 
that  has  been  thoroughly  researched  and  documented  during  the 
past  10-12  years  (1,9).  This  research  project  will  attempt 
to  establish  a  model  for  space  diversity  gain  by  using  well 
known  linear  diversity  combining  techniques  and  giving  con¬ 
sideration  to  the  spatial  separation  of  the  earth  stations. 

Standards 

The  criterion  used  to  measure  the  performance  of  each 
signal  selection  scheme  will  be  the  probability  of  bit  error 
(P^)  versus  the  average  signal-to-noise  ratio  (SNR)  per  bit, 
Yb,  of  the  received  downlink  signal.  This  measurement  is 


the  standard  criterion  for  comparing  the  performances  of 
digital  communications  systems.  These  values,  (?e  vs.  y^) * 
will  be  calculated  and  plotted  for  each  selection  scheme 
during  the  computer  simulation.  These  results  will  allow 
for  a  direct  comparison  of  all  of  the  developed  signal  selec¬ 
tion  schemes  and  a  measurement  of  the  benefits  derived  from 
using  space  diversity  transmission  schemes  as  compared  to 
the  result  from  the  use  of  a  standard  satellite  system  with¬ 
out  space  diversity. 

Approach 

A  detailed  investigation  of  possible  linear  diversity 
combining  techniques  for  optimum  signal  selection  will  be 
conducted  to  determine  which  methods  are  best  suited  for  this 
satellite  communications  application.  For  each  diversity 
technique,  the  criteria  for  determining  the  optimum  received 
signal  will  be  fully  developed.  After  a  complete  develop¬ 
ment  and  analysis  of  each  of  the  proposed  diversity  signal 
selection  methods,  a  computer  simulation  will  model  the  oper¬ 
ational  environment  and  provide  a  performance  measurement 
for  each  diversity  technique. 


II .  Satellite  Background  Information 

The  major  elements  in  a  satellite  communications  system 
include  the  satellite,  a  network  of  earth  stations,  and  multi 
pie-access  communications  equipment  by  which  many  earth  sta¬ 
tions  can  operate  through  a  single  satellite  (22*131).  A 
satellite  communications  link  consists  of  a  satellite  orbit¬ 
ing  in  space  and  one  or  more  earth  stations.  The  performance 
of  this  communications  link  is  greatly  influenced  by  the 
earth-space  propagation  channel.  One  of  the  major  differ¬ 
ences  between  satellite  links  and  terrestrial  telecommunica¬ 
tions  is  the  distances  that  the  carriers  must  be  transmitted. 
Since  the  strength  of  the  radiated  signal  diminishes  with  the 
square  of  the  distance  it  travels,  the  satellite  signals  are 
relatively  weak  when  received  at  the  distant  end.  However, 
the  amount  of  received  carrier  power  determines  the  ability 
of  the  receiver  to  demodulate  or  decode  the  information 
(11*83).  In  satellite  systems  it  is  extremely  important  to 
know  the  key  parameters  that  directly  determine  this  received 
power  so  that  these  systems  can  be  properly  analyzed.  Thus, 
an  understanding  of  the  basic  power  flow  equations  associated 
with  satellite  channels  is  necessary. 

Antenna  Gain 

The  antennas  used  in  satellite  links  do  not  radiate 


equally  in  all  directions  but  are  designed  to  focus  the  radi¬ 
ation  in  a  particular  direction.  The  increase  in  power 


achieved  by  focusing  the  antenna  is  referred  to  as  the  gain 
of  the  antenna.  An  isotropic  antenna  which  radiates  equally 
in  all  directions  has  an  antenna  gain  =  1 .  The  gain  of  the 
antenna,  Gt(0z,0^)  a  functi°n  of  the  angular  direction  of 
the  receiver,  where  (0  ,0.)  refers  to  the  azimuth  and  eleva- 

^  Xj 

tion  angle,  respectively,  measured  from  a  coordinate  system 
centered  at  the  transmitting  antenna  (11:83).  Therefore, 
antennas  can  be  characterized  by  their  gain  pattern  G.  (0_,0„) 
which  indicates  how  the  antenna  gain  is  distributed  spatially 
with  respect  to  the  antenna  coordinate  system. 

The  most  important  parameters  of  an  antenna  pattern  are 
its  gain  (the  maximum  value  of  the  gain  pattern) ,  its  beam- 
width  (a  measure  of  the  angle  in  which  most  of  the  gain 
occurs),  and  its  sidelobes  (amount  of  gain  in  the  off-axis 
direction)  (11:87).  For  most  satellite  communications  pur¬ 
poses,  antenna  patterns  should  be  highly  directional  with 
most  of  the  gain  concentrated  over  a  narrow  beamwidth  and  a 
very  small  part  of  the  radiated  energy  in  the  sidelobes. 
However,  the  gain  and  beamwidth  is  dependent  on  the  type  of 
antenna  used . 

The  gain  of  a  transmitting  antenna  is 

Power  received  from  the  antenna 


Power  received  if  the  radiation  were  isotropic 


Although  the  antenna  gain  is  dependent  on  the  type  of 
antenna  used,  a  useful  approximation  for  this  gain  is 


where 


=  aperture  cross-sectional  area  of  transmitting 
antenna 

Z  =  the  efficiency  of  the  antenna  aperture 
X  =  wavelength  of  transmission 

X,  the  wavelength  of  the  radiation,  is  inversely  proportional 
to  its  frequency  f* 

X  =  ~y~  (2.3) 

where  C  is  the  velocity  of  light.  Thus,  the  antenna  gain 
approximation  can  be  written  such  that 

Z4rrf  2A 

G  =  - = - -  (2.4) 

CT 

where 

C  =  2.99  x  10^  meters/second 
It  is  obvious  that  the  gain  of  the  antenna  is  dependent 
on  the  frequency  of  transmission  and  the  efficiency  of  the 
antenna  aperture . 

EIRP 

The  power  of  the  transmitter  is  often  referred  to  in  a 
manner  which  takes  the  gain  of  the  antenna  into  consideration. 
It  is  referred  to  as  the  effective  isotropic  radiated  power 
(EIRP) .  EIRP  is  the  power  of  a  transmitter  and  isotropic 
antenna  that  would  achieve  the  same  result  as  the  transmitter 
and  antenna  in  question. 


where 


Pfc  =  available  antenna  input  carrier  power  from 
L  the  transmitter  power  amplifier,  including 
circuit  coupling  losses  and  antenna  radia¬ 
tion  losses 

Gt{*z,<V  =  gain  of  transmitting  antenna 

Now  we  can  define  the  received  carrier  power  collected 
by  the  receive  antenna  that  is  normal  to  the  transmitter  as 


P 


r 


(EIRP)L  A 
_ _  a  r 

4ttZ2 


(2.6) 


where 


L  =  atmospheric  loss 

d 

A  =  the  aperture  of  the  receive  antenna 

=  UZ/4rT)Gr(0z’,5y )  (2.7) 

where 


Gr(0z’,0  ’)  =  azimuth  and  elevation  angles  of  the  trans¬ 
mitter  relative  to  the  receiver  coordinate 
system 

X  =  carrier  wavelength 


An  equivalent  expression  for  Eq  (2.6)  is 

Pr  =  (ElRP)L_LpGr(0z’ ,0£')  (2.8) 

where  we  define 

Lp  =  (  X/4ttD)  2  =  (C/4iTfD)2 

Lp  is  called  the  propagation  or  free  space  loss  that  is 
incurred  by  the  communications  link.  It  is  dependent  on  the 
frequency,  f,  and  the  distance,  D. 
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Eq  (2.7)  can  also  be  restated  in  a  more  useful  form 


where 

(P  ) dB  =  (EIRP)dB  +  (L  )dB  +  (L_) dB  +  (G  JdB  (2.9) 

r  p  a  r 

External  Noise 

In  addition  to  receiving  the  desired  signal,  the  re¬ 
ceiving  antenna  also  gathers  other  forms  of  electromagnetic 
energy  that  may  be  present.  These  other  radiations  may  be 
interference  or  noise  that  can  distort  the  transmitted  sig¬ 
nal.  The  major  external  noises  to  a  satellite  communications 
link  are: 

1 .  Atmospheric  Losses/Noise 

2 .  Thermal  Noise 

3 .  Sky  Noise 

4.  Galactic  Noise 

5.  Cosmic  Noise 

The  combination  of  all  sources  of  interference  form  a 
total  noise  level  for  the  receiving  system  of  the  satellite 
link  (16:123)  . 

Atmospheric  Losses 

There  are  many  frequency  dependent  effects  on  the  propa¬ 
gation  channel,  in  addition  to  the  free  space  loss,  that 
increases  the  total  transmission  loss  of  the  communications 
link  (6:162).  These  propagation  effects  can  be  categorized 
as  atmospheric  losses.  The  primary  causes  of  atmospheric 
losses  for  satellite  communications  links  are: 


1 .  uncondensed  water  vapor 

2 .  molecular  oxygen 

3.  rainfall 

Water  Vapor  and  Molecular  Oxygen  Attenuation.  These 
two  losses  are  relatively  constant  and  easily  calculated, 
although  the  attenuation  levels  vary  within  the  carrier  fre¬ 
quency.  The  main  reason  for  the  attenuation  due  to  water 
vapor  and  molecular  oxygen  is  that  some  of  their  associated 
atmospheric  particles  have  molecular  resonances  in  the  micro- 
wave  frequency  band. 

Below  10  GHz,  atmospheric  losses  are  minimal,  (usually 
less  than  2  dB) ;  however,  at  higher  frequencies,  the  attenua¬ 
tion  begins  to  increase  rapidly.  The  higher  levels  of  atten¬ 
uation  occur  at  the  frequencies  having  wavelengths  corre¬ 
sponding  to  specific  gas  molecules  in  the  atmosphere  (11s 94). 
Therefore,  the  attenuation  caused  by  water  vapor  and  oxygen 
varies  greatly  with  the  frequency.  Severe  attenuation  due  to 
water  vapor  molecules  occurs  at  ~22  GHz  while  the  severe 
attenuation  due  to  oxygen  molecules  occurs  at  ~60  GH.  Thus, 
the  frequencies  for  the  RF  carrier  for  satellite  links  are 
chosen  in  the  frequency  range  where  these  types  of  attenua¬ 
tion  are  minimized. 


Rainfall  Attenuation  Effects.  The  effects  of  rainfall 
attenuation  on  satellite  communications  varies  with  frequency 
however,  it  can  have  a  severe  effect  on  signal  strength  at 
frequencies  above  10  GHz.  Rainfall  is  the  most  serious  atmos 
pheric  effect  on  a  satellite  link;  therefore,  special 


consideration  must  be  given  to  the  effects  of  precipitation 
on  the  earth-space  propagation  path  when  designing  a  satel¬ 
lite  communications  system  (8*2-1).  At  frequencies  above 
10  GHz,  rain  is  known  to  cause  significant  attenuation, 
absorption,  scattering,  amplitude  and  phase  scintillations, 
depolarization,  and  bandwidth  decoherency.  All  of  these 
phenomena  combine  to  degrade  the  reliability  of  the  satel¬ 
lite  communications  link.  Its  obvious  impact  on  the  satel¬ 
lite  system  is  cost,  since  methods  must  be  incorporated  into 
the  system  to  combat  rainfall  effects  if  uninterrupted 
service  is  desired  (10:191). 

Rainfall  is  a  complicated  subject  primarily  because  of 
the  complexity  of  the  rainfall  itself.  Rainfall  can  assume 
many  different  forms  depending  on  the  geographical  location, 
time  of  year,  time  of  day,  temperature,  altitude,  and  other 
earth  station  site  characteristics  (8:2-1). 

Types  of  Rain.  As  described  in  Ref  8,  chapter  2,  there 
are  three  basic  types  of  rain:  stratiform,  convective,  and 
cyclonic.  In  mid-latitude  regions  of  the  world,  stratiform 
rain  is  a  type  of  rain  that  primarily  occurs  in  the  spring 
and  fall.  It  exhibits  stratified  horizontal  extents  of  hun¬ 
dreds  of  kilometers  (km) ,  vertical  heights  of  4-6  km,  with 
durations  exceeding  one  hour  and  rain  rates  less  than  about 
25  mm/hr  (1  inch/hr).  Stratiform  rain  has  a  rain  rate  and 
duration  which  often  necessitates  the  use  of  link  power  mar¬ 
gin  to  exceed  the  attenuation  caused  by  a  25  mm/hr  rain  rate 


Convective  rains  are  a  result  of  vertical  atmospheric 
motions  producing  vertical  transport  and  mixing.  The  con¬ 
vective  flow  occurs  in  a  cell  whose  horizontal  extent  is 
usually  several  kilometers.  The  cell  usually  extends  to 
heights  exceeding  the  freezing  layer  because  of  convective 
upwelling.  The  cell  may  be  isolated  or  embedded  in  a  thunder 
storm  associated  with  a  weather  front.  Because  of  the  motion 
of  the  front  and  the  sliding  motion  of  the  cell  along  the 
front,  the  duration  of  high  rain  rate  is  usually  only  several 
minutes .  This  type  of  rain  is  the  most  common  source  of  high 
rain  rates  in  the  U.S.  and  Canada. 

Tropical  cyclonic  storms  (hurricanes)  occasionally  pass 
along  the  U.S.  eastern  and  southern  coasts  during  August  to 
October.  These  circular  storms  are  typically  50-200  km  in 
diameter,  move  at  10-20  km/hr,  extend  upward  to  8  km,  and 
have  high  rain  rates  that  exceed  25  mm/hr. 

Stratiform  and  cyclonic  rainfall  cover  large  geographic 
areas  and  the  spatial  distribution  of  total  rainfall  from 
such  a  storm  is  assumed  to  be  uniform.  Therefore,  the  rain 
rate  averaged  over  several  hours  will  be  similar  for  earth 
stations  located  tens  of  kilometers  apart  (8:2-2).  Convec¬ 
tive  storms  are  localized  and  generally  produce  non-uniform 
distributions  of  rainfall  and  rain  rate  for  a  given  storm. 

The  total  rainfall  and  rain  rate  varies  significantly  over  a 
distance  of  10  km  (8:2-2). 

There  are  several  proposed  models  for  modeling  the  rain¬ 
fall  attenuation  effects  on  satellite  communications. 


However,  all  of  these  models  contain  frequency  dependent 
coefficients  and  require  the  rain  rate  statistics  for  the 
particular  earth  station  sites  since  the  rain  rates  vary  with 
location.  Currently,  rain  rate  statistics  are  developed  by 
data  collection  experiments  at  the  desired  location. 

Receiver  Noise 

In  addition  to  the  external  noise  collected  by  the 
antenna,  noise  is  also  generated  within  the  components  of  the 
receiver  immediately  following  the  antenna.  The  contribu¬ 
tions  from  these  sources  of  interference  combine  to  define  a 
total  noise  level  for  the  receiving  system.  This  noise  level 
set  the  minimal  required  power  from  the  desire  transmitter 
to  achieve  reliable  communications. 

Noise  Temperature.  The  amount  of  receiver  noise  present 

.  O 

is  defined  by  the  receiver  noise  temperature,  T  e^.  This 
parameter  is  an  effective  equivalent  temperature  that  an 
external  noise  source  would  need  to  produce  to  have  the  same 
amount  of  receiver  noise  (11:101). 

T°  =  T°  +  ( F- 1 )  290°  (2.10) 

eq  b 

where 

T°,  =  background  noise  temperature  accounting  for  con- 

tj  tributions  collected  by  the  antenna  from  external 
sources 

F  =  noise  figure  of  receiver 

A  noise  source  of  temperature  T°eg  produces  an  effective 
noise  spectrum  level  at  the  antenna  input  of 


Watts/hz 


(2.11) 


where 

-23 

K  =  Boltzmann's  constant  =  1.379  x  10  W/fc-Hz 

=  -228.6  dBW/k-Hz 

The  total  noise  entering  the  receiver  over  a  bandwidth, 

bff-  1s 


where 


B 


RF 


P  = 
n 


KT°eqBRF 


bandwidth  of  the  receiver 


(2.12) 


Therefore,  the  total  receiver  noise  can  be  calculated  from 
the  knowledge  of  the  T0^  of  the  receiver. 

Noise  Figure.  The  internal  thermal  noise  is  accounted 
for  in  T°  by  specifying  the  receiver  noise  figure,  F,  which 
depends  on  the  specific  electronic  circuit  components  follow¬ 
ing  the  antenna.  The  noise  figure  for  an  electronic  system 
is  defined  as 


Total  system  output  noise  power 

F  =  -  (2.13) 

Output  system  noise  power  due  to 
system  input  noise  at  290°K 

When  defined  as  this,  ( F  —  1 )  290°K  is  the  temperature  of  an 
equivalent  noise  source  at  the  system  input  that  produces 
the  same  contribution  to  the  system  output  noise  as  the 
internal  noise  of  the  system  itself  (11:103).  For  an  example 
on  the  specifics  on  calculating  F,  see  Ziemer  and  Tranter 
(23:478-81) . 


One  of  the  most  useful  indicators  of  the  performance  of 
a  communications  receiver  is  the  ratio  of  the  received  car¬ 
rier  power  to  the  total  noise  power  of  the  receiver.  This 
RF  carrier-to-noise  ratio  (CNR)  is  defined  as: 


(2.14) 


where 


Pr  -  <EIRP)LpLaGr 
Pn  *  KT°eqBeq 

The  CNR  indicates  the  relative  strength  of  the  desired  trans¬ 
mitter-receiver  power  and  the  total  noise  interference 
(11:104).  Eq  (2.14)  shows  how  the  different  link  parameters 

affect  receiver  CNR.  Since  the  receiver  bandwidth  is 

RF 

often  dependent  on  the  modulation  scheme,  it  is  convenient 
to  isolate  the  RF  link  power  parameters  by  normalizing  the 
CNR  with  respect  to  bandwidth.  The  resulting  RF  carrier-to- 
noise  level  ratio  is  defined  as 


(C/Nq)  =  <EIRp)LpLaGr  (2.15) 

Obviously  (c/Nq)  is  now  independent  of  .  In  a  digital  com¬ 
munications  system,  the  (C/NQ)  ratio  allows  us  to  compute 
directly  the  receiver  bit  energy-to-noise  level  ratio  as 

Eb/Nc  =  (C/N0)Tb  (2.16) 

where 

Tfc  =  bit  time 


By  knowing  the  link  (C/NQ) ,  which  depends  only  on  the  RF 
link  parameters,  we  can  compute  the  analog  CNR  by  dividing 
by  Brf  or  the  digital  Eb/NQ  by  multiplying  by  the  bit  time. 

It  is  convenient  to  express  Eg  (2.15)  as 

(C/N0)  =  (EIRP/k)  (LpLa)  (Gj/r°eq)  (2.17) 

where 

(EIRP/K)  =  transmitter  parameter 
(LpLa)  =  propagation  parameter 

(Gr/T°eg)  =  receiver  parameter 

Thus,  each  subsystem  of  the  link  contributes  to  the 
overall  (C/Nq) .  The  only  effect  of  the  receiving  system  is 
the  ratio  Gr/T°eq.  Therefore,  as  far  as  the  receiver  is  con¬ 
cerned,  performance  can  be  maintained  with  a  lower  receiver 

gain  (smaller  antenna)  if  T°  (noise)  can  be  reduced. 

eg 

There  is  a  direct  trade-off  of  receiver  antenna  size  and 
receiver  noise  temperature  in  achieving  a  desired  performance 

Power  Link  Equations 

It  is  convenient  to  distinguish  between  the  uplink  and 
downlink  transmission  paths  and  their  corresponding  RF  car- 
rier-to-noise  level  ratios.  The  uplink  eguation  can  be  des¬ 
cribed  as 

(C/N  )  =  EIRP)  -L  +  G  -T-K-L  dB-Hz  (2.18) 

o'u  '  e  pu  r  au 

where 
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EIRP)  =  earth  station  EIRP,  dBW 
e 

Lpu  =  free  space  path  loss,  uplink,  dB 

Gr  =  gain  of  the  satellite  receive  antenna,  dBi 

T  =  satellite  receiver  noise  temperature,  dB-K 
K  =  Boltzmann’s  constant  =  -228.6  dBW /K -Hz 
L  =  uplink  atmospheric  losses,  dB 

all 

The  downlink  equation  is  described  as 


(C/N0)u  -  £IRP)S  -  BOo  -  V  +  Gr  -  T  -  K  -  Lad  <2'19> 


dB-HZ 


where 


EIRP)  =  satellite  station  EIRP,  dBW 


BO  =  output  backoff,  dB 


L  ,  =  free  space  path  loss,  downlink,  dB 


pd 

G 


r 

T 

K 


=  gain  of  the  earth  station  receive  antenna,  dBi 


Jad 


earth  station  receiver  noise  temperature,  dB-K 
Boltzmann’s  constant  =  -228.6  dBW/K-Hz 
downlink  atmospheric  losses,  dB 


Note:  The  output  backoff  term  is  expressed  in  dB  from  the 
saturated  output  point  relative  to  saturated  power.  Since 
the  power  transfer  function  is  nonlinear  for  the  satellite 
transponder,  BO  is  a  nonlinear  function  (2:14). 


Digital  Signaling  Methods  for  Satellite  Communications  Systems 
In  modern  satellite  communications  systems,  the  trend 


has  been  towards  digital  instead  of  analog  communications . 

A  digital  communication  system  transmits  information  wave¬ 
forms  that  are  first  converted  into  a  sequence  of  data  sym¬ 
bols  which  are  then  transmitted  over  the  satellite  link  by 
encoding  onto  RF  carriers,  (11:34),  as  depicted  in  Figure 
2.1.  A  major  advantage  of  digital  communications  is  that 
information  can  be  transmitted  relatively  error-free  with 
less  carrier  power  than  is  required  for  an  analog  system. 

For  satellite  systems,  digital  communications  has  the  addi¬ 
tional  advantage  that  transmitted  carrier  waveforms  can  be 
carefully  controlled  in  terms  of  amplitudes  and  frequency 
spectra,  simplifying  satellite  hardware  design  (11:34). 

The  digital  information  to  be  transmitted  over  the  chan¬ 
nel  is  assumed  to  be  a  sequence  of  data  symbols  occurring  at 
a  uniform  rate  of  R  symbols/second.  For  binary  signaling, 
the  set  of  data  symbols  is  restricted  to  a  "1”  bit  and  a 
"0"  bit  and  each  symbol  may  be  transmitted  directly  by  send¬ 
ing  a  waveform  corresponding  to  a  "l"  or  a  waveform  corres¬ 
ponding  to  a  "0". 

As  an  alternative  to  transmitting  the  binary  bits  di¬ 
rectly,  the  information  sequence  can  be  grouped  into  blocks, 
where  each  block  consists  of  k  bits  (19:140).  With  k  bits/ 

V  t 

block,  there  are  2  =  M  distinct  blocks.  Therefore,  M  dif¬ 

ferent  waveforms  are  needed  to  transmit  the  k-bit  blocks 
unambiguously . 


Encoded 

Baseband 

Waveform 


The  commonly  used  measure  of  performance  for  digital 
communications  systems  is  the  relationship  between  bit  error 
rate  and  signal-to-noise  ratio  (SNR)  (2t38).  The  SNR  for 
digital  systems  is  usually  expressed  as  the  ratio  energy  per 
bit  to  the  noise  spectral  density  (EbAJQ) .  Generally,  the 
desired  modulation  technique  will  be  the  one  that  satisfies 
the  required  system  constraints  and  requires  the  least  value 
for  Ej_)/Nq  for  a  specified  threshold  error  rate. 

Digitals  systems  can  be  easily  compared  on  the  basis  of 
ideal  bit  error  rate.  In  comparing  digital  satellite  com¬ 
munications  systems  the  ideal  bit  error  rate  is  used  to  deter¬ 
mine  the  best  possible  performance.  The  following  assumptions 
are  made  for  the  "ideal  bit  error  rate"  conditions  (2:39): 

1.  No  phase  or  frequency  uncertainties  are  introduced 
by  earth  station  or  satellite  oscillators. 

2.  No  forward  error  correcting  codes  are  used. 

3.  Carrier  and  clock  recovery  is  perfect,  so  that 
there  is  no  recovered  carrier  phase  or  clock  jitter. 

4.  The  only  source  of  errors  is  the  addition  of  random 
additive  white  gaussian  noise  by  the  channel. 

Multiphase  Signaling  Waveforms.  In  a  general  M-ary 
signaling  system,  the  possible  transmitted  signal  waveforms 
may  be  denoted  as  lsm(t)(,  m  =  1, 2,  ...,  M.  We  shall  consider 
these  waveforms  to  be  bandpass  signals,  and  represent  them  as 


s  (t)  =  Re[u  (t)exp(jnf  t)] 


(2.20) 


where 


m  =  1,2,  .  .  .M 

i  um ( t ) }  =  the  equivalent  low-pass  waveforms 
Re  =  the  real  part  of 

The  M  signals  are  characterized  individually  by  their  energy, 
defined  as, 

T 

i  =  S  S  2(t)  dt 
c  m  0  m 

T 

=  0.5  S  I  a  (t)l 2  dt  (2.21) 

o  m 

where 

m  —  1,2,... M 

Note:  We  shall  assume  that  all  M  signaling  waveforms  have 

equal  energy  such  that 

£  =  £  (2.22) 

m 

where 

m  1,2,  ..  .M 

For  M-ary  phase  shift  keying  (PSK),  the  baseband  sig¬ 
naling  waveforms  }sm(t)|  are  used  to  phase  modulate  the  RF 
carrier.  The  general  representation  for  a  set  of  M-ary  phase 
signaling  waveforms  is 

sm(t)  =  Re{u(t)  exp  j  ( 2TTf  ct  +  (  2tt  (m-1 )  Al) +X)(  (2.23) 

where 

m  —  1,2,  ...,  M 
0  <  t  <  T 

X  is  a  fixed  initial  phase 


u(t)  =  0 


otherwise 


then  sm(t)  may  be  expressed  as  a  PSK  signal 

s  (t)  =  Acos[2-rTf  t  +  (2rr  (m-1) /M) +X]  (2.24) 

m  c 

where 

m  —  1  ,  2  ,  «  .  *  M 
0  <  t  <  T 

Performance  of  M-ary  PSK.  The  probability  of  bit/symbol 
error  analysis  is  thoroughly  presented  in  Proakis,  (Ref  19, 
Chap  4) .  The  results  of  this  analysis  ares 

Case  Is  M  =  2,  PSK  Waveforms  =  Antipodal  Waveforms 

pb =  q£<2V0,5]  (2,25) 

where 

Pfe  =  probability  of  bit  error 
Yb  =  Eb/N0  =  SNR  Per  bit 

Q  (x)  =  1  -  /^(I/Utt)  °*5)  exp(-B2/2)  dB  (2. 25. a) 

Case  2s  M  =  4,  X  =  tt/4 

Pm  =  2Q[(2Yb)0,5]  -  [Q(2Yb)0,5]2  (2.26) 

where 

P  =  probability  of  k-bit  symbol  error 


(2.27) 
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Case  3:  For  M  >  8  and  y  >  10  dB 


Pm  =  2Q((2KYb)0,5  sin  (-tt/M)  ) 


P,  =  P_/k 
b  m 


(2.28) 


Differentially  Encoded  PSK  (DPSK) .  Ordinary  PSK  re¬ 
quires  coherent  processing  and  accuracy  constraints  on 
estimating  the  phase  increase  as  M  increases.  A  procedure 
which  allows  for  easier  phase  and  symbol  synchronization  is 
called  differential  PSK.  In  M-ary  DPSK,  instead  of  trans¬ 
mitting  signals  with  phases  corresponding  to  the  M  possible 
k-tuples,  we  transmit  phases  corresponding  to  the  change 
(difference)  between  the  previous  k -tuple  and  the  current 
k-tuple  (4:V-65) . 

Performance  of  DPSK. 


For  M  =  2: 


Pfa  =  0.5  exp ( - Y b ) 


For  M  =  4: 


Pfe  =  Q (a, b) -0 .5IQ( ab) exp  [0.5(a2  +  b2)] 


(2.29) 


(2.30) 


where 


7  7  K 

Q(a,b)  =  exp[0.5(a  +  b  )  ]]  T.  (a/b)  I  (ab) ,  for  b  >  a  >  0 

K  =  0  K 


Ia(x)  =  I 


IxZ2) 


a  +  2k 


,  X  >  0 


K=0  K!  r(a  +  K  +  1) 


•--.v 
•*  /V*  1 
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III.  Rainfall  Characteristics 


Rainfall  has  a  significant  effect  on  satellite  communica¬ 
tions  at  frequencies  above  10  GHz.  Due  to  the  scattering 
effect  by  rain  which  causes  attenuation,  depolarization,  and 
interference,  the  characteristics  of  rain  became  an  impor¬ 
tant  parameter  to  satellite  communications  systems  designers. 
In  order  to  understand  and  appreciate  the  problem  of  rainfall 
attenuation,  it  is  vital  to  discuss  the  basic  characteristics 
of  rain  and  its  associated  statistics . 

Rainfall  Physical  Structure  Characteristics 

The  variability  of  rain  with  respect  to  space  and  time 
has  been  documented  as  early  as  the  1950’s  (1:369).  Measure¬ 
ments  on  the  ground  have  shown  that  differences  in  rainfall 
intensities  can  occur  within  distances  of  a  few  hundred  feet. 
The  understanding  of  the  formation  and  distribution  of  large 
raindrops  and  high  rainfall  rates  is  of  great  importance  to 
the  designers  of  communications  satellite  systems. 

Heavy  rain  usually  occurs  with  cellular  structures  of 
limited  horizontal  and  vertical  extent  (1:369).  An  important 
characteristic  of  raincells  is  their  tendency  to  form  into 

groups  of  cells  within  a  region  called  a  small  mesoscale  area 
.  2 

which  can  be  100-400  km  .  As  many  as  seven  cells  have  been 
observed  in  a  small  mesoscale  area  and  up  to  six  small  meso- 

scales  areas  within  a  large  mesoscale  area  which  is  about  1000 

2  .  ... 
km  .  The  variations  of  the  cell  distribution  with  a 
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mesoscale  structure  are  not  purely  random;  the  form  of  the 
mesoscale  pattern  can  be  greatly  influenced  by  local  climato¬ 
logical,  meteorological,  or  topographical  effects  (1j369). 

The  lifetime  of  rain  showers  and  thunderstorms  is  longer 
than  those  of  an  individual  cell;  however,  the  lifetime  of  a 
single  cell  can  vary  greatly  from  as  little  as  a  few  minutes 
to  as  long  as  several  hours.  The  dynamics  of  raincells 
appear  to  fall  into  two  basic  categories:  those  in  which 
the  rainfall  intensity  within  the  cell  is  well  preserved  over 
periods  of  minutes  while  the  cell  itself  moves  across  the 
terrain,  and  those  in  which  the  cells  are  basically  station¬ 
ary  but  exhibit  large  fluctuations  in  rainfall  intensity  with 
time  ( 1  s  369)  . 

Rain  Rates 

Rain  is  normally  measured  as  the  rate  (mm/hr)  at  which 
the  rainfall  is  collected  by  a  horizontal  surface  near  the 
ground  at  the  location  of  interest.  Therefore,  in  order  to 
obtain  the  rain  rate  of  a  particular  location,  an  experiment 
designed  to  obtain  these  statistics  must  be  performed  at  that 
location.  Since  rain  varies  greatly  with  time,  the  raw  data 
contains  too  much  information  to  be  useful,  thus  the  data  is 
processed  to  reveal  the  long-term  statistical  characteristics 
of  rain.  The  most  important  of  these  is  the  cumulative  dis¬ 
tribution  function  (CDF)  (10:199),  which  shows  the  fraction 
of  time  a  certain  rain  rate  is  not  exceeded  (see  Figure  3.1). 
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NUMBER  or  MINUTES  per  year  rain  rate  exceeds  abscissa 


In  addition,  for  each  rain  rate,  a  typical  distribution 
of  drop  sizes  in  a  volume  of  space  above  and  near  the  rain 
gauge  can  be  determined  as  a  function  of  the  observed  instan¬ 
taneous  rain  rate  (10:193).  Although  several  rainfall  dis¬ 
tributions  have  been  developed,  some  of  the  more  common  aver¬ 
age  dropsize  distributions  used  in  rain  attenuation  calcula¬ 
tions  are  (19:321): 

1.  Laws  and  Parsons  Distribution:  This  tabular  dis¬ 
tribution  has  been  found  to  be  a  reasonable  measure  for  an 
average  dropsize  spectrum  in  continental  temperate  rainfall. 
It  is  probably  the  most  widely  used  distribution  currently 
available . 

2.  Marshall-Palmer  Distribution:  This  negative  exponen¬ 
tial  distribution  is  a  fairly  good  fit  for  the  average  drop- 
size  spectrum.  It  is  most  applicable  to  the  widespread  rain 
in  continental  temperate  climates  although  it  has  a  tendency 
to  overestimate  the  number  of  small  drops. 

3.  "Thunderstorm"  Distribution  of  Joss  et.  al.t  This 
negative  exponential  distribution  was  fitted  to  the  average 
dropsize  spectrum  measured  in  convective  rain. 

4.  "Drizzle"  Distribution  of  Joss  et.  al . :  This  nega¬ 
tive  exponential  distribution  was  obtained  by  fitting  the 
average  dropsize  spectrum  of  very  light  widespread  rain  com¬ 
posed  mostly  of  small  drops. 

A  typical  Laws  and  Parson  dropsize  distribution  of  rain¬ 
drop  diameters  and  velocities  for  various  rain  rates  is  shown 
in  Table  3.1. 


Measured  Raindrop  Terminal  Velocities  and  Dropsize  Distributions  at  Ground  Leve] 
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Rain  Scattering 

In  order  to  calculate  a  specific  attenuation  factor 
that  corresponds  to  a  specific  rain  rate,  one  must  understand 
the  phenomena  of  rain  scattering.  When  rain  particles  enter 
the  path  of  the  satellite  communications  downlink  signal,  the 
particles  both  absorb  and  scatter  the  radiation  associated 
with  the  transmitted  signal.  The  interaction  between  the 
rain  particles  and  the  electromagnetic  field  associated  with 
the  signal  can  be  satisfactorily  computed  by  using  electro¬ 
magnetic  scattering  theory.  However,  exact  solutions  for  the 
scattering  problem  are  only  available  for  simple  shapes  and 
distributions  of  the  dielectric  properties  of  the  scatterer 
within  the  scatterer  (7:176).  Rain  attenuation  may  be  rea¬ 
sonably  modeled  using  spheres  with  the  assumptions  that  water 
is  a  homogeneous  lossy  dielectric,  the  drops  are  randomly 
distributed  with  a  volume  in  accordance  with  a  Poisson  proc¬ 
ess,  and  the  drop  size  distribution  is  known  (6:460).  More 
complex  models  of  raindrops  have  been  developed  to  account 
for  the  nonspherical  physical  description  of  rain  but  these 
models  have  only  been  partially  successful  and  only  at  lower 
rain  rates  (7:176). 


Rain  Attenuation  Models 

The  most  fundamental  quantity  in  the  calculation  of  rain 
attenuation  statistics  for  space-earth  paths  is  the  attenua¬ 
tion  per  unit  distance,  A.  There  are  two  general  approaches 
commonly  used  to  calculate  A,  a  theoretical  technique  and  an 


empirical  method  based  upon  the  approximate  relationship 
between  A  and  the  rain  rate  R.  Both  methods  rely  on  the 
solution  of  the  scattering  problem  by  Mie  (17:377,  7:176, 
18:318-9).  The  validity  of  the  theoretical  approach  is  well 
established  (18:318)  and  its  calculations  have  been  conducted 
for  a  wide  variety  of  frequencies.  However,  the  complexity 
of  this  approach  makes  it  unsuitable  for  use  by  satellite 
communications  systems  designers. 

The  theoretical  analysis  of  the  rain  attenuation  results 
in  the  expression 

OO 

A  =  4.343  /  Q  (D)  N(D)dD  (3.1) 

o  t 

where 

Q t  =  (4rrA2)  Re[s(0,D)] 

S  (0,D)  =  forward  scattering  amplitude 

N(D)dD  =  number  density  of  raindrops  with  equivalent 
diameter  D  in  the  interval  dD 

Although  it  was  believed  that  this  equation  included  only  the 
effects  of  single  scattering,  it  also  included  the  effects  of 
all  the  forward  multiple  scattering  processes  (i.e.,  those 
processes  which  all  paths  between  scatterers  have  components 
in  the  forward  direction)  (18:319).  As  suggested  by  Hogg  and 
Chu  (14:1310-2),  these  are  the  only  multiple  scattering  proc¬ 
esses  which  significantly  contribute  to  the  coherent  rain 
attenuation . 

The  second  method  of  calculating  the  rain  attenuation  is 
based  on  the  approximation 
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where 


a,b  =  functions  of  frequency,  raindrop  size  distribution, 
and  rain  temperature 

This  expression  serves  as  a  good  approximation  for  the  full 
Mie  solution  obtained  in  the  theoretical  approach  for  fre¬ 
quencies  greater  than  1  GHz  and  is  suitable  for  most  design 
applications.  The  difficulty  with  this  method  is  obtaining 
accurate  values  for  the  coefficients  a  and  b.  Many  methods 
for  obtaining  the  coefficients  have  been  developed  including 
linear  (b=l)  and  nonlinear  relations,  values  for  the  coef¬ 
ficients  at  different  frequencies  by  fitting  the  results  of 
Mie  calculations,  and  experimentally  obtained  empirical 
values  for  a  and  b  (18:318).  As  expected,  there  are  minor 
discrepancies  in  the  values  for  these  coefficients.  An  often 
used  set  of  coefficient  values  is  presented  by  Olsen,  Rogers, 
and  Hodges  (18).  This  method  presents  a  complete  set  of 
values  of  a  and  b  for  frequencies  1-1000  GHz  where  the  values 
of  specific  attenuation  obtained  at  several  rain  rates  were 
used  to  determine  a  and  b  by  means  of  logrithmic  regression 
for  different  combinations  of  frequency,  dropsize  distribu¬ 
tion,  and  temperature  (See  Tables  3.2,  3.3,  3.4). 

Statistical  Behavior  of  Rain 

In  addition  to  the  models  to  calculate  a  rain  attenua¬ 
tion  factor,  there  has  been  attempts  to  characterize  the 
behavior  of  rain  attenuation.  One  such  effort  is  the  work  by 


TABLE  3.2 

Regression  Calculations  For  a  and  b  in 
A  =  aR*5  (dBAm)  as  Functions  of  Frequency 
and  Dropsize  Distribution  with  Rain  Temperature 


=  20°C 


r - —  r 

l\ 

HP 

J-T 

J-D 

3-04.10"? 

3.17.10*? 

5.H.I0'? 

2.12.10*? 

5.17.10'? 

1.52*10*: 

6.75*10*; 

1.14*10*! 

5  53*10  : 

1 . 14.10  ■ 

I.SJ.IOlJ 

1.15.10'; 

2.02*10*! 

7-ic.io*; 

2.03.10  , 

2- *2. 10  * 

1 . 73.10*i 

3-19*1 0*! 

8.29*10*; 

3.20.10*! 

3. 57.10*: 

2-39.10*; 

4  04.10*! 

1.01.10*! 

4.67*10*1 

4.99*10  ! 

3.11*10*; 

‘  *0*10*! 

1.01.10*! 

6.46*10*! 

4.75.10  , 

3.7S.1C  , 

i.4fe.io*; 

2.45.10*; 

8.59*10*! 

1.12*10  : 

5.15*10*! 

1  38.10*: 

7.49.10*! 

1.40.10'; 

1.29*10  , 

1.06*10*: 

2  13*10* J 

2.17*10’: 

2.13.10'; 

2 .82*10  , 

2.04*10  i 

3. 26*10*: 

4.59*10*: 

3.04.10'; 

4.34*10  : 

3.78.10*: 

4.84.10*; 

1.01.10*; 

4.25*10*: 

5.73*10*: 

9.44.10*: 

6.74.10*; 

7  10*10'; 

1.55*10*; 

9.20*10*, 

l.ii.io*; 

1,01.10*; 

lido’; 

7.52.10'i 

1.27.10*; 

1.67.10*; 

1.37*10*; 

2.42.10*; 

1.47.10'; 

i.9B*io  ; 

7.33*10** 

i,0i*)o*; 

3.25.10'; 

1.22*10'; 

328.10’; 

4. *9.10*; 

8  59.10** 

3-57.10'; 

5.04.10'* 

2.20.10, 

|.83«10’Z 

7.51*10** 

0.103 

4.56*10  ; 

0.113 

0.143 

0.127 

0.101 

7.76*10*Z 

0  160 

0.220 

0.191 

0.273 

0.118 

0.235 

0.337 

0.269 

0.360 

0.167 

0.312 

0.452 

0.  360 

0.439 

0  225 

0. 404 

0.648 

0.572 

0.610 

0.368 

0.652 

0  775 

0.796 

0.771 

0.541 

0  802 

0  650 

1.01 

0.014 

0.735 

0  034 

0.902 

i.ii 

0.810 

0.936 

1  03 

0.938 

i.»» 

0  876 

1.14 

1.09 

0  958 

1  53 

0  995 

1.33 

1.13 

0  972 

1.64 

1.09 

1.50 

1.17 

0.982 

1.74 

111 

1.66 

1.33 

1.02 

1.96 

t -07 

2  05 

1.50 

1.06 

2.09 

0  993 

2  47 

1  46 

l  07 

2.17 

0.991 

2.69 

1.41 

0  993 

2.74 

0  997 

2.84 

1.40 

0.980 

2  22 

0  981 

2  9? 

1  36 

0  958 

7  18 

0  959 

2.95 

1  31 

0  976 

7  14 

0.937 

2  97 

1  27 

0  902 

2  11 

0  917 

2  94 

1.25 

0  006 

7  08 

0  897 

2  94 

1.23 

0  *71 

2  05 

0  081 

2  44 

>  r 

0  461 

2  03 

0  067 

7  41 

1  1  20 

0.851 

2  01 

C.856 

7  «  1 

l  ■  ,i 


l 


Source:  (18:323) 


Vi.%  ' j  *  *  •  »  *  .  •  »  •  .  •  . 


For  a  and  b  in 
>ns  of  Frequency 
tin  Temperature 


1 

B 

J-T 

*+ 

on? 

0.9*8 

0.15* 

0.901 

0.843 

0.909 

0.978 

0.171 

0.970 

0.847 

0.931 

1.01* 

0.892 

1.005 

0.C&7 

0.95$ 

1.053 

0-917 

1.190 

8.858 

0.991 

1.094 

0.944 

1.2*3 

0.894 

1.008 

1.134 

0.971 

1 .171 

8.870 

1.033 

i.m 

0.998 

1.270 

8.877 

1.077 

1.214 

1.0*1 

1.77* 

8.881 

1.105 

1.222 

1.068 

l.ltO 

0.WI 

1.119 

1.21* 

1.083 

1.145 

O.t'O 

i.m 

i.m 

1.090 

1.111 

0.*1* 

i  m 

1.185 

1.097 

1.099 

8.928 

i.m 

i.m 

1.092 

1.0W 

0.937 

1.124 

1.161 

1.091 

1.077 

8.938 

1.111 

1.1S2 

1.090 

1.064 

0.942 

1.119 

1.135 

1.088 

1.077 

8.954 

1.111 

1.096 

1.082 

0.951 

0.171 

1.091 

1.037 

1.065 

0.887 

0.988 

1  059 

0-969 

1.037 

0.826 

1.004 

1  071 

0.910 

1.007 

0.784 

1.01* 

0.911 

0.867 

0.968 

0.761 

1.025 

0.910 

0.818 

0.90* 

0.714 

1.077 

i  0.859 

0.796 

0.854 

0.964 

1.016 

0.817 

0  785 

0  817 

0.659 

0.997 

0.787 

C  781 

0  789 

0.673 

0.974 

0.785 

0.777 

0.764 

0.668 

0.953 

0.74? 

0  774 

0.741 

0.646 

0.933 

0.739 

0.771 

0.774 

0  674 

0.914 

0  733 

0.770 

0.717 

0.611 

0.899 

0711 

0.769 

C.6A6 

0.610 

0.846 

0  874 

0.765 

0.65? 

0.609 

0.754 

0.859 

0.763 

0  638 

0.608 

0.739 

0.860 

0.764 

0  618 

0  61? 

0.704 

0.658 

0  764 

0  61? 

0.610 

0.684 

0  658 

0  765 

0.6H 

0.612 

0.969 

0.864 

0.767 

0.611 

0  611 

0.85S 

0  668 

0  768 

0  608 

0.609 

0.834 

0.671 

0  769 

0.609 

0.609 

0.622 

0.67? 

0.770 

-0610 

0.611 

0.617 

0.673 

0  770 

0610 

0  813 

0.814 

9.674 

0.771 

0-911 

0.9U 

0  613 

Source:  (18:324) 


Lin  (15 t 557-81)  regarding  the  behavior  of  rain  attenuation. 

Lin  gathered  31  sets  of  experimental  data  on  the  statistics 
of  microwave  rain  attenuation  at  frequencies  above  10  GHz  in 
the  U.S.,  England,  Japan,  Italy,  and  Canada  and  proposed  many 
ideas  dealing  with  the  statistical  behavior  of  rain  attenua¬ 
tion. 

The  summary  of  his  findings  are  as  follows: 

1.  the  distribution  of  rain  attenuation,  a(t),  in  dB 
is  approximately  lognormal  within  the  attenuation  range  of 
1-50  dB. 

2.  the  distribution  of  the  fade  duration,  f(A),  is  also 
approximately  lognormal . 

Lin's  theory  on  the  lognormal  behavior  of  rain  attenua¬ 
tion  suggests  that  three  parameters  which  depend  on  geographic 
locations  determine  the  distribution  of  the  rain  attenuation. 
These  three  parameters  are: 

P  =  probability  of  rainfall,  the  expected  fraction  of 
time  that  rain  falls  at  the  location  of  the  satel¬ 
lite  earth  station 

a  =  standard  deviation  of  the  attenuation,  a(t) 
a 

a  =  median  of  the  attenuation,  a(t) 
m 

Lin  also  stated  that  am  is  frequency  dependent  and  increases 
almost  linearly  with  the  path  length  because  the  median  rain 
is  usually  small  and  almost  uniform  over  the  entire  path.  c?a 
decreases  slightly  as  frequency  increases  and  decreases 
slightly  as  the  path  length  increases  because  of  the  aver¬ 
aging  effect  of  the  propagation  volume. 
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The  equation  describing  the  lognormal  distribution  of 
rain  attenuation  is 

P(a>A)  =  0.5(Po)  (0.5+Q[(log10A-tiin)/(20,5aa)])  (3.3) 

where 

PQ  =  probability  of  rainfall 

=  standard  deviation  of  log^Q  a  during  the 
a  raining  time 

“a  "  lo9l0an. 

=  mean  value  of  log^a  during  the  raining  time 

Similarly,  the  distribution  of  the  fade,  f(A),  is 
described  by 

P  (f (a)>b  )  =  0.5(0.5+Q[(logl0b-uf)/20*5Of)])  (3.4) 

where 

(i ^  =  mean  of  log1Q 

=  standard  deviation  of  log^CffA)] 

Due  to  the  significance  of  rain  attenuation  to  the  satel¬ 
lite  communications  systems  designer,  an  ideal  method  that 
would  permit  an  accurate  extrapolation  of  measured  rain  char¬ 
acteristics  to  different  geographical  locations  or  other  time 
periods  is  needed.  Since  detailed  measurements  of  the  various 
important  statistics  of  rainfall  exist  for  only  a  few  locations 
throughout  the  world,  it  would  be  of  great  value  if  these 
statistics  could  be  accurately  applied  or  converted  to  other 
areas,  especially  those  with  different  climates.  However,  a 


40 


good  deterministic  theory  of  rain  generation  and  rainfall 
has  not  been  developed,  therefore  it  is  not  possible  to 
apply  it  to  the  prediction  of  localized  rain  occurrences . 


I 


Diversity  Techniques  For  Satellite  Communications 


With  the  recent  growth  of  satellite  communications, 
allocations  in  the  higher  frequency  bands  have  been  required 
to  relieve  the  congestion  in  the  8/7  GHz  band.  One  such 
allocation  is  the  30/20  GHz  band,  which  is  being  considered 
for  many  government  applications  (9:1-1).  Unfortunately,  a 
severe  problem  with  rain  attenuation  exists  in  the  30/20  GHz 
frequency  band  (as  noted  in  Sections  II  and  III) . 

There  are  many  techniques  that  can  be  employed  to  over¬ 
come  the  rain  attenuation.  In  order  to  overcome  the  deep 
fades  on  the  uplink  or  downlink,  satellite  transmitters  could 
employ  higher  power  levels  or  the  associated  earth  station 
could  use  a  larger  antenna  (1:369).  However,  another  concept 
used  to  combat  the  effects  of  rain  attenuation  is  diversity, 
which  is  employed  to  reduce  the  margin  required  for  rain 
attenuation.  The  main  concept  of  diversity  is  that  the  prob¬ 
ability  of  an  error  can  be  reduced  by  using  two  or  more 
received  signals  to  determine  the  actual  transmitted  signal 
instead  of  one. 

Description  of  Basic  Diversity  System 

There  are  three  basic  types  of  diversity:  angle,  fre¬ 
quency,  and  space  (or  site) .  The  use  of  angle  diversity 
requires  two  or  more  satellites  which  are  separated  physi¬ 
cally  such  that  the  fading  on  any  earth-space  path  has  a  low 
correlation  with  the  fading  on  other  paths  (9:3-1) 


(See  Figure  4.1).  Due  to  the  requirement  for  two  satellites, 
angle  diversity  is  very  costly  and  it  has  been  shown  to  be 
not  as  effective  as  the  other  two  diversity  methods  (9*1—2)  . 
Frequency  diversity  requires  that  the  satellite  be  equipped 
with  two  transponders,  one  operating  in  a  high  frequency  band 
and  the  other  in  a  lower  band.  When  the  rain  attenuation  is 
minimal  the  higher  band  is  used  and  the  lower  frequency  band 
transponder  is  used  when  the  rain  attenuation  becomes  sig¬ 
nificant.  Although  the  fades  in  the  two  bands  are  correlated, 
the  fade  depths  in  the  lower  frequency  band  are  less  than  the 
higher  frequency  band.  Frequency  diversity  incurs  a  relativ¬ 
ely  high  cost  for  the  additional  satellite  transponder  but 
the  performance  has  been  rated  the  highest  of  the  three  diver¬ 
sity  methods  (9:1-2). 

However,  the  focus  of  this  thesis  effort  is  on  space 
diversity  techniques  that  can  be  used  to  combat  rain  attenua¬ 
tion  in  the  downlink  of  a  satellite  communications  link. 

The  use  of  site  diversity  in  the  downlink  requires  that  two 
or  more  earth  terminals  be  physically  separated  such  that  the 
fading  on  any  space-earth  path  has  a  low  correlation  with  the 
fading  on  other  paths  (9:2-1).  Dual  diversity  is  the  most 
basic  type  of  space  diversity  but  higher  orders  of  diversity 
may  be  needed  at  locations  with  a  history  of  heavy  rains. 
However,  the  higher  the  order  of  diversity,  the  more  costly 
the  system  due  to  the  additional  earth  stations. 

A  block  diagram  of  a  dual  downlink  space  diversity  system 
is  shown  in  Figure  4.2.  Since  the  diversity  equipment 
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combines  the  digital  outputs  of  the  two  earth  stations,  the 
configuration  is  independent  of  whether  Frequency  Division 
Multiple  Access  (FDMA)  or  Time  Division  Multiple  Access  (TDMA) 
is  used  on  the  satellite  communications  link  (9:2-1).  The 
diversity  equipment  could  be  collocated  with  one  of  the 
earth  stations  or  it  could  be  located  at  a  separate  site. 

The  optimum  spacing  between  the  earth  stations  is 
dependent  on  the  characteristics  of  the  rain  in  that  partic¬ 
ular  area.  In  practice,  optimum  spacing  may  not  be  achiev¬ 
able  because  of  limitations  on  the  availability  of  sites  for 
earth  stations.  However,  the  spacing  will  usually  be  in 
excess  of  one  mile  necessitating  the  use  of  a  terrestrial 
link  to  connect  at  least  one  earth  station  to  the  diversity 
equipment . 

Bit  Stream  Synchronization.  When  using  diversity  on  the 
downlink,  there  will  usually  be  a  difference  in  the  transmis¬ 
sion  path  lengths  from  the  satellite  to  the  combiner.  How¬ 
ever,  due  to  the  satellite  motion,  this  difference  in  path 
lengths  will  not  remain  constant.  The  path  length  difference 
and  its  variation  are  functions  of  the  relative  geographical 
locations  of  the  earth  stations  and  the  satellite  (9:2-3). 
Therefore,  due  to  the  path  length  difference,  some  adjustment 
or  control  must  be  made  to  maintain  bit  synchronization  of 
the  bit  streams  into  the  combiner.  If  the  path  length  dif¬ 
ference  is  a  significant  fraction  of  a  bit,  a  common  prac¬ 
tice  is  to  use  a  first-in-first-out  (FIFO)  buffer  on  the 
shorter  path  as  shown  in  Figure  4.2  so  that  the  bits  arrive 
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at  the  combiner  in  synchronism  (9s 2-5)  .  If  the  variation  in 
the  path  length  difference  is  a  significant  fraction  of  a 
bit,  a  timing  tracking  loop  is  needed  to  maintain  synchro¬ 
nism.  In  many  cases,  this  loop  operates  similarly  to  the 
code  tracking  loop  in  a  pseudonoise  receiver  except  that  one 
of  the  two  bit  streams  is  used  as  a  reference  (9s2-5). 

Linear  Diversity  Combining  Techniques.  When  using 
spatial  diversity  a  combining  technique  must  be  used  to 
choose  the  optimum  signal  from  the  diversity  branches.  The 
most  common  methods  for  combining  diversity  channels  outputs 
ares 

1.  Selection  combining.  Selection  of  the  channel  with 
the  highest  signal  level.  Measurement  of  the  signal  level 
is  required  and  this  method  can  use  either  hard  or  soft 
decisions . 

2.  Maximum  Selection  combining.  Selection  of  the 
received  symbol  with  the  largest  amplitude.  The  method 
requires  soft-decision  outputs. 

3.  Equal  Gain.  Summing  the  soft-decision  outputs  of 
each  channel . 

4.  Square  Law  combining.  Squaring  the  soft-decision 
output  of  each  channel  followed  by  linear  summing. 

5.  Maximum  Ratio  combining.  Weighting  the  soft-deci¬ 
sion  output  of  each  channel  by  the  SNR.  The  measurement  of 
the  SNR  is  required. 

The  diversity  improvement  is  dependent  on  many  factors 
including  the  rain  fading  characteristics  of  that  particular 
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area  (1:369).  However,  theoretically,  maximum  ratio  com¬ 
bining  is  the  most  efficient,  and  selection  combining  is  the 
least  ( 9 s 2  —  8 ) .  Since  the  difference  in  performance  between 
the  different  combining  methods  is  relatively  small,  the 
decision  of  which  method  to  use  can  vary  with  the  particular 
implementation  depending  on  the  desired  performance  criteria 
and  cost  restrictions. 

Space  Diversity  Measurements 

The  rain  attenuation  effects  on  the  transmitted  signal 
for  a  slant-path  link  can  be  measured  directly,  or  inferred 
from  the  measurement  of  other  related  parameters  by  using 
one  of  the  following  techniques  (1:369-70): 

1.  direct  reception  of  a  beacon  signal  transmitted 
from  a  geostationary  satellite  at  the  required  frequency. 

2.  radiometric  measurements  using  either  the  sun  as  a 
source  or  the  apparent  change  in  sky  temperature  with  rain¬ 
fall  intensity. 

3.  rapid-response  raingauges . 

4.  radar. 


As  noted  in  Table  4.1,  there  are  advantages  and  limita¬ 
tions  for  each  of  the  different  measurement  techniques  for 
obtaining  slant-path  attenuation.  The  three  major  advantages 
of  the  radiometer  are  relative  cheapness,  extreme  reliability 
and  straightforward  data  reduction.  For  gathering  long-term 
data  for  slant-path  attenuation,  the  most  cost-effective 
solution  is  the  radiometer  (1:370).  Fast-response  raingauges 


Measurement  Techniques  For  Obtaining  Slant— Path  Attenuation 


although  inexpensive,  do  not  provide  data  that  can  be  easily- 
transformed  into  slant-path  attenuation  data.  In  order  to 
provide  an  accurate  transformation  of  collected  data  into 
slant-path  attenuation  data,  an  extensive  array  of  rain- 
gauges,  accurate  telemetry  and  large  computing  facilities 
are  required,  which  greatly  increases  the  associated  costs 
of  this  measurement  technique.  Radar  is  a  very  useful  tool 
for  investigating  propagation  phenomena;  however,  there  are 
also  some  potential  problem  areas.  Unless  great  care  is  used 
in  eliminating  all  echoes  due  to  solid  or  solid/liquid  phase 
particles,  very  large  errors  will  accrue  in  the  calculation 
of  the  slant-path  attenuation  (1:370).  The  best  measurement 
method  for  obtaining  the  desired  slant-path  propagation  data 
is  the  satellite  beacon  technique.  However,  to  orbit  radio 
beacons  covering  a  vast  range  of  frequencies  is  extremely 
expensive  thus  the  opportunity  of  conducting  satellite-beacon 
propagation  experiments  rarely  exists  except  for  a  limited 
period  at  a  few  frequencies  and  locations  (1:370). 


Diversity  Performance  Assessment  Techniques 


There  are  two  major  methods  currently  employed  in  assess¬ 
ing  the  performance  of  space  diversity  systems:  diversity 
gain  and  diversity  improvement.  Diversity  gain  is  defined  as 
the  difference  between  the  path  attenuation  associated  with 
a  single  earth  station  and  the  diversity  modes  of  operation 
for  a  given  percentage  time  (13:280).  It  uses  (percentage  of) 
the  time  as  a  parameter  and  expresses  the  diversity 


improvement  as  the  difference  in  decibels  which  is  the  effec' 
tive  increase  in  fade  margin  between  the  single  earth  sta¬ 
tion  site  and  the  diversity  system  fading  statistics.  Diver 
sity  improvement  is  defined  as  the  ratio  of  the  percentage 
time  a  given  fade  depth  is  exceeded  on  a  single  path  to  the 
percentage  time  that  the  same  fade  depth  is  exceeded  for  the 
diversity  system  (13:280).  An  example  of  these  two  perform- 


observation  is  probably  due  to  the  fact  that  the  calculation 
of  the  improvement  factor  involves  division  by  a  percentage 
time  which  is  often  very  small  and  therefore  subject  to  con¬ 
siderable  uncertainty  (13*281).  Since  diversity  gain  is  a 
function  of  the  single  earth  station  site  fading  depth  char¬ 
acteristics,  it  is  commonly  used  to  compare  sets  of  space- 
diversity  measurements  conducted  over  different  periods  of 
time  and  during  different  weather  conditions.  Hence,  the 
diversity  gain  method  of  performance  measurement  for  diver¬ 
sity  systems  is  generally  preferred  to  the  diversity  improve¬ 
ment  method . 

Much  research  has  been  conducted  to  develop  a  model 
which  accurately  predicts  the  performance  of  a  diversity  sys¬ 
tem  with  respect  only  to  the  separation  distances  between  the 
earth  stations  and  not  to  the  carrier  frequency  (1,13).  One 
of  the  best  known  models  of  this  type  is  the  empirical  rela¬ 
tionship  for  path  diversity  gain  developed  by  D.B.  Hodge  in 
the  mid-70's  (13*280-1).  Hodge  gathered  diversity  gain  data 
for  millimeter  wavelength  earth-space  propagation  paths  in 
Ohio  and  New  Jersey.  Although  these  measurements  were  col¬ 
lected  in  different  locations  and  over  different  time  periods, 
the  resulting  diversity  gain  data  were  very  consistent. 

Hodge  used  this  to  develop  an  empirical  relationship  for 
diversity  gain  as  a  function  of  separation  distance  and  sin¬ 
gle  earth  station  fade  depth.  He  recognized  that  the  diver¬ 
sity  gain  data  followed  the  same  relationship  for  each 


value  of  the  single  earth  station  fade  depth.  This  rela¬ 
tionship  for  diversity  gain,  Gd,  is 

Gd  =  a(l  -  exp(-bD))  (4.1) 

where 

a  =  A  -  3.6(1  -  exp(-0.24A)) 
b  =  0.46(1  -  exp (-0.26 A)) 

A  =  single  earth  station  fade  depth  in  dB 
D  =  earth  station  separation  distance  in  km 

The  coefficients  a  and  b  depend  only  upon  the  single  earth 
station  fade  depth  and  were  found  by  performing  a  minimum  rms 
error  fit  for  each  fade  depth. 

Although  Hodge's  relationship  for  diversity  gain  accu¬ 
rately  modeled  the  data  measurements  collected  for  his  experi¬ 
ment,  there  have  been  other  similar  experiments  where  the 
collected  data  were  not  accurately  modeled  by  the  Hodge  rela¬ 
tion.  Fimbel  and  Juy  assumed  that  diversity  gain  was  fre¬ 
quency  dependent  and  recalculated  the  coefficients  a  and  b 
from  Hodge's  formula  and  were  able  to  develop  a  relationship 
that  more  accurately  predicted  the  results  of  the  data  meas¬ 
urements  that  they  had  collected  in  the  United  Kingdom  (1:374). 

In  an  effort  to  improve  its  accuracy,  Hodge  later  modi¬ 
fied  his  model  for  diversity  gain  as  a  function  of  earth 
terminal  separation  distance,  link  frequency,  elevation 
angle,  and  angle  between  the  baseline  and  the  path  azimuth 
(12:1393-9).  His  new  model  for  diversity  gain,  GD,  is 


where 

Gd  =  function  defined  by  Eq  (4.1) 

Gf  =  function  of  frequency 

Ge  =  function  of  elevation  angle 

G^  =  function  of  orientation  of  the  baseline  relative 
to  the  propagation  path 

Gf  =  1.64  exp ( -0 .025f )  (4.3) 

where  f  =  link  frequency  in  gigahertz 

Ge  =  0.00492e  +  0.834  (4.4) 

where  e  =  elevation  angle 

Gb  =  0.00177b  +  0.887  (4.5) 

where  b  =  orientation  of  the  baseline  relative  to  the 
propagation  path  in  degrees 

A  set  of  thirty-four  experiments  were  conducted  to  val¬ 
idate  the  Hodge's  improved  diversity  model  and  the  model 
reproduced  the  original  data  set  with  a  rms  error  of  0.73  dB. 

The  research  in  this  area  has  continued  but  so  far  no 
single  model  developed  has  been  able  to  accurately  predict 
all  of  the  diversity  gain  data  measurements  gathered  by 
experiments  conducted  throughout  the  world.  It  is  widely 
believed  that  the  difficulty  in  establishing  a  general  model 
for  predicting  the  diversity  gain  with  respect  to  separation 


distance  is  that  there  are  several  causes  of  variation  for 
the  diversity  gain  data  measurements.  The  three  basic  cate¬ 
gories  which  cause  diversity  gain  variation  are  climatic, 
geographic,  and  systemic  (1:372). 

The  climatic  factors  which  cause  variations  in  diver¬ 
sity  gain  are  associated  with  large  scale  weather  patterns. 
They  are  rain  lifetimes,  rain  cell  sizes,  rain  cell  separa¬ 
tion,  and  rain  cell  orientation.  In  different  locations  or 
during  different  periods,  these  climatic  factors  may  vary 
and  cause  a  variation  in  the  diversity  gain  measurements  with 
respect  to  earth  station  separation  distances. 

The  geographic  causes  relate  to  the  fact  that  there  can 
be  local  variations  in  climate  within  a  limited  area  as 
small  as  10  km.  The  geographic  feature  that  seem  to  enhance 
or  reduce  the  rainfall-intensity  probability  at  given  loca¬ 
tions  is  generally  referred  to  as  the  microclimate  of  a 
region  (1:373).  These  local  variations  in  climate  are  of 
particular  importance  for  space  diversity  communications  sys¬ 
tems  since  knowledge  of  where  the  high  and  low  probability 
areas  of  intense  rainfall  would  greatly  aid  in  the  siting 
of  earth  stations. 

The  systemic  causes  of  variations  in  diversity  gain 
measurements  are  variation  with  frequency,  variation  with 
elevation  angle,  variation  with  baseline  orientation,  and 
variation  with  site  separation.  The  first  three  variations 
deal  with  the  concept  that  rain  attenuation  varies  with 
frequency  and  the  slant- path  of  the  rain  with  respect  to  the 


earth  station.  The  fourth  variation,  separation  distance, 
deals  with  the  fact  that  a  good  site  location  is  dependent 
on  several  factors,  including  the  microclimate,  prevailing 
wind  direction,  and  the  viewing  bearing.  Thus,  the  separa¬ 
tion  distance  is  dictated  by  the  results  of  these  other  fac¬ 
tors.  However,  site  separations  of  approximately  10  Km 
seem  to  provide  diversity  gain  close  to  the  optimum  although 
some  regions  will  require  larger  separation  distances  to 
achieve  the  same  diversity  gain  because  of  geographic  and 
climatic  reasons  (1:374). 
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V.  Analysis  of  Selected  Linear  Diversity  Combining  Techniques 

In  the  previous  section,  it  was  mentioned  that  a  diver¬ 
sity  combining  technique  must  be  implemented  within  the 
diversity  network  in  order  to  successfully  achieve  spatial 
diversity  improvements.  In  this  section,  some  of  the  diver¬ 
sity  combining  techniques  will  be  examined  so  that  a  proba¬ 
bility  of  bit  error  analysis  can  be  performed  on  the  selected 
combinations  of  diversity  techniques  and  digital  signaling 
schemes.  The  diversity  techniques  to  be  examined  are  the 
selection  diversity,  equal  gain,  and  maximal  ratio  techniques 
and  the  digital  signal-schemes  are  restricted  to  DQPSK  and 
B-ary  PSK.  Therefore,  the  possible  combinations  of  these 
diversity  techniques  and  signaling  schemes  will  result  in  six 
different  probability  of  bit  error  expressions. 

Considerations  for  Bit  Error  Analysis 

In  performing  such  an  analysis,  there  are  assumptions 
that  must  be  stated  in  order  to  validate  the  analysis.  Some 
of  the  assumptions  for  the  development  of  the  statistics  for 
the  two  digital  signaling  schemes  are: 

1.  The  channel  through  which  the  signaling  waveforms 
are  transmitted  is  assumed  to  have  no  bandwidth  restrictions 
(i.e.,  the  ideal  case). 

2.  The  signals  transmitted  through  the  channel  are  cor¬ 


rupted  by  additive  white  gaussian  noise  (AWGN) . 


3.  The  receivers  at  the  earth  stations  are  synchro¬ 


nized  to  the  satellite  transmitter. 

4.  The  carrier  phase  is  Known  perfectly  at  the  receiver 
in  order  to  facilitate  ideal  coherent  demodulation. 

5.  All  of  the  lovpass  signal  representations  in  each 
particular  signaling  scheme  have  the  same  energy  and  are 
equiprobable . 

The  diversity /signaling  scheme  combinations  will  be 
compared  on  the  basis  of  average  bit  error  rate,  Pg.  If 

Pg ( Yb^  =  conditional  probability  of  bit  error  given  y^ 
for  a  specified  signaling  scheme 

PD(Yb^  =  probability  of  density  function  of  yfc  for  a 
particular  diversity  combining  technique 

then  [as  in  Eq  7.4.14  from  (Ref  19:474)J 

PB  *  -C  VVb>PD(Yb)d*b  (5'1) 

Consider  first  DQPSK,  Referring  to  Section  II  Eq  (2.30) 
of  this  thesis 

PB,DQPSK(Yb)  =  Q (a,b)  -  0.5IQ(a,b)  exp  [-0.5(a2  +  b2)] 
where 

a  =  1.0824  (0.5yb)°’5 
b  =  2.6131  (0.5yb)0'5 

Q(a,b)  =  exp(_0*5(a2  +  b2))I  (|)K  I„(ab)f  for  b  >  a  >  0 

K- 0  b  K 
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A  further  simplification  is  provided  by  Shanmugam 


(21:425).  For  4-ary  DPSK 


PB,DQPSK^Yb^  ~  2Qt(4Ybsin  (  /g) )  ]  (5.2) 


For  M-ary  PSK  signaling  the  probability  of  a  symbol 
error  is  expressed  by 

TT 

PM  =  1-  f  M  .p(e)d0  (5.3) 

-rr 
M 

where 


p(0)  =  i  exp  ( -yb)  [  1  +  ( 4TTYb)  °  ‘ ' 5  ) 

(2tt) 


(2Vbcose)0-5 

•exp(ybcos20)  S  exp ( -0 . 5X2 ) dx] 

— oo 


The  probability  of  a  bit  error  for  an  8-ary  PSK  system 
with  Gray  code  assignment  of  3-tuples  (Figure  5.1)  can  be 
expressed  by 

VsPSK^b*  '  l  ,(K  errors)  (peK>  (5-4 

K  —  1 


where 

PeK  =  probability  of  error  for  K  bit  errors 
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1,1 


Figure  5.1  Gray  Code  for  8-ary  PSK 


Therefore, 


PB,  8PSK 


3tt/8  5tt/8  7tt/8 

(YoJ  =  C  S  p(6)d6  +  2  J’  p(0)d0  +  S 

u  tt/8  3tt/8  5tt/8 


9tt/8  llrr/8 

+2  S  p(0)d6  +3  S  p( 0)  d0  +  2  J 
7 tt/8  9rr/8  1 


15tt/8 

+  S  p(0)d0 

13tt/8 


Since  p(0)  is  an  even  function,  Eq  (5.5)  reduces  to 

3  tt/8  5tt/8  7tt/ 

Pr  rpsk(yJ  =  2  S  p(0)d0  +  4  S  p(0)d0  +  4  S 
B,  8PSK  b  <rr/8  3tt/8  5tt/8 


tt 

+  /  p(0)d0 


Due  to  the  complexity  of  evaluating  the  substitution 


of  Eq  (5.3)  into  Eq  (5.6),  the  approximation  denoted  by 
Eq  (2.28)  is  commonly  used  for  the  probability  of  bit  error 
for  high  order  M-ary  PSK  systems.  For  8-ary  PSK 


PB,8PSK^Yb^  K  3  sin(rr/8)3 


(5.7) 


Assumptions  for  Diversity  Analysis 

Having  developed  the  probability  of  bit  error  for  each 
signaling  scheme,  we  must  now  look  at  the  different  linear 
diversity  techniques.  Since  the  digital  signaling  techniques 
being  used  require  coherency,  we  will  limit  the  diversity 
techniques  to  those  which  are  most  commonly  used  with  coher¬ 
ent  signaling  schemes.  These  three  diversity  techniques, 
selection,  equal  gain,  and  maximal  ratio,  can  all  be  gener¬ 
ally  described  by  the  equation 


f(t)  =  a1f1(t)  +  ...  +  aLfL(t)  =  a^fj(t)  (5.8) 

where 

f(t)  =  selected  output  of  diversity  system 

a.  =  combining  coefficient  proportional  to  the  gain 
J  of  the  jth  channel 

fj  =  received  signal  of  the  jth  channel 

L  =  number  of  diversity  channels 


There  are  several  assumptions  about  the  channel  char¬ 
acteristics  that  must  be  stated  before  attempting  further 
development.  The  basic  assumptions  are  made  for  the  channels 
when  employing  the  selected  linear  diversity  techniques  ares 


1.  The  channel  impulse  response  will  be  modeled  as  a 
zero  mean  complex-valued  gauss ian  process  so  that  its  enve¬ 
lope  at  any  time,  t,  is  Rayleigh-distributed. 

2.  The  amplitude  variations  in  the  received  signal, 

i.e.,  signal  fading,  are  caused  by  the  time-variant  multi- 

path  channel  characteristics.  Due  to  assumption  1,  the 

channel  attenuation,  a,  is  Rayleigh-distributed;  therefore 
2 

a  has  a  chi-square  probability  distribution  with  2  degrees 
.  .  2 

of  freedom.  Since  is  dependent  on  a  ,  is  also  chi- 
square-distributed . 

3.  The  information-bearing  signal  transmitted  through 
the  channel  is  large  in  comparison  to  the  bandwidth  of  the 
transmitted  signal  so  that  the  channel  is  considered  fre- 
quency-nonselective . 

4.  The  signal  bandwidth  and  the  signaling  interval 
shall  be  chosen  small  enough  so  that  the  channel  attenua¬ 
tion  and  phase  shift  are  essentially  fixed  for  the  duration 
of  at  least  one  signaling  interval;  therefore  the  channel 
can  be  considered  slowly  fading. 

5.  The  channel  fading  is  sufficiently  slow  so  that  the 
phase  shift  can  be  perfectly  estimated  from  the  received 
signal . 

6.  The  noise  in  each  channel  is  independent  of  the 
signal  and  additive;  and  the  received  signal,  fj(t)  =  Sj(t) 

+  nj(t),  where  s^  and  n^  are  the  signal  and  noise  components, 
respectively,  in  the  jth  channel. 
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TV  '  •  *.J  V  V  ■  i 


7.  The  noise  components, 
and  have  zero  means . 


n  j  ( t) 


are  locally  un correlated 


Probability  of  Error  in  Rayleigh  Channel  with  No  Diversity 
In  order  to  evaluate  the  benefits  of  a  diversity  sys¬ 
tem,  the  probability  of  a  bit  error  must  be  calculated  for 
DQPSK  and  8-ary  PSK  in  a  f requency-nonselective,  slowly 
Rayleigh  fading  channel  without  diversity  using  the  same 
assumptions  about  the  channels  as  for  the  diversity  case. 

By  assumption  2,  is  chi-square  distributed,  therefore, 

pno  diversity**!,1  =  ^  exp  ‘"Vb^b1  <5'9) 

vb 

Using  Eq  (5.1),  the  probability  of  a  bit  error  can  be  obtained 
for  DQPSK  and  8-ary  PSK  in  a  Rayleigh  channel  without  diver¬ 
sity  . 

For  DQPSK, 

00  1 

Pb(DQPSK,  no  diversity)  =  f  —  exp  (~Yb/vb) 

°  Yb 


'PB, DQPSK (Yb)dYb 

where 

PB,DQPSK^Yb^  =  2Q^4vbsin  (tt/8)  )  ] 

For  8-ary  PSK, 

00  _ 

PB(8PSK,  no  diversity)  =  /  —  *exp(-Yb/Yb) 

°  Yb 


(5.10) 


’PB,8PSK(Yb^dYb 


(5.11) 
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where 


PB,8PSK^Yb^  ”  3  ®^6Yb)  sin(TT/8)3 
Selection  Diversity 

The  basic  principle  of  the  selection  diversity  tech¬ 
nique  is  that  at  any  given  time,  the  system  picks  out  the 
best  of  the  L  noisy  signals  f  ,  f 2,  ,  fL,  and  uses  that 

one  alone;  and  the  other  signals  do  not  contribute  to  f(t) 
(3sl081).  The  ideal  selection  diversity  combiner  is  defined 
as  choosing  for  the  system  output,  during  each  "instant"  of 
time,  the  signal  from  that  receiver  which  has  the  largest 
SNR  (See  Figure  5.2) .  Therefore,  the  parameter  used  to 
determine  the  optimum  diversity  channel  is  received  SNR. 

If  K  denotes  the  index  of  a  channel  for  which  E,  ,,  >  E.  ., 

bK  bj' 

j  =  1,2,  ...  L;  where  j  =  K,  then  this  type  of  system  can 
be  characterized  by  the  Eq  (5.8)  where 

1,  for  j  =  K 

a.  =  (5.12) 

J  0,  for  j  =  K 

Several  methods  can  be  employed  for  selecting  the  optimum 
diversity  channel  signal.  One  of  the  easiest  is  if  all  the 
receivers  have  the  same  noise  level  then  the  channel  with 
highest  total  instantaneous  power  is  selected. 

Analysis  of  Selection  Diversity 

Given  the  assumption  of  short-term  fading,  the  perform¬ 
ance  of  selection  diversity  can  be  obtained  from  the 


Figure  5.2  Selection  Diversity 

statistics  of  the  SNR  at  the  selector  output.  In  partic¬ 
ular,  the  p.d.f.  of  y^,  the  instantaneous  SNR  in  the  Kth 
channel  is  needed.  Following  the  development  of  Schwartz, 
et  al.,  (20:434-5),  can  be  easily  obtained.  Assum¬ 

ing  Rayleigh  fading  and  defining  y  as  the  average  of  y, 
over  the  short-term  fading  for  the  Kth  channel  so  that 

,  ,  eK2I  g 

*<*  =  EhcKl  =  '  -I.  -  (5-13) 

I\ 

where 

aK  =  attenuation  factor  of  Kth  channel 
$  =  energy  of  signal  transmitted 
Nk  =  power  spectral  density  of  Kth  channel 
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yields  the  individual  p.d.f.'s 


PD^bK1  =  exP<-VbK/VcK>  <5'14) 

cK 

and  the  associated  cumulative  distributions  are 

P(YbK<x)  =  l-exp(-X/ycK)  (5.15) 

Over  short-term  fading,  all  the  events  in  which  the 
selector  output  SNR,  Yb>  is  less  than  X  is  exactly  the  set 
of  events  in  which  each  of  the  channel  SNR's  is  simultane¬ 
ously  below  X  (20:435) .  Since  the  fading  is  assumed  inde¬ 
pendent  in  each  of  the  L  channels,  the  probability  that 
Yb  <  X  is  the  product  of  the  independent  probabilities  that 
each  ybK  <  X. 

Therefore, 

L  L 

P(Yb<X)  =  n  P(YbK<X)  =  n  [ 1 -exp ( -X/y cK ) 3  (5.16) 

K~  1  K~  1 

The  corresponding  p.d.f.  for  is 

pD<V  ■  aMi  :i-^p(-x/ycK>3 } x  -  vb 

n  ~  1 

=  df;{K51i:i-exp<-vb/vcK)]} 

If  we  assume  all  diversity  channels  have  equal  mean  SNR  over 
the  short-term  fading,  ycK  =  Yc  »  then 


(5.17) 


(5.18) 
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P(vb<x)  =  [l-exp(-X/Yc)]i 


(5.19) 


pD<Yb)  =  o4"  exP  (-Yb/vc)[l-exp(-vb/?c)]L'1  (5.20) 


Now  the  probability  of  bit  error  can  be  calculated  for 
the  selection  diversity  combining  technique  using  DQPSK  and 
8-ary  PSK . 

For  selection  diversity  using  DQPSK, 


P„ (DQPSK ,  selection  diversity)  =  £  —  exp  (-y./y  ) 

B  o  Yc  D  c 


[l-exp(-Yb/vc)] 


'PB, DQPSK (Yfa)dYb(5.2l) 


where 


PB,DQPSK^Yb^  2Q[(4Ybsin  (tt/8)  )  ] 


For  selection  diversity  using  8-ary  PSK, 


Pr(8PSK,  selection  diversity)  =  J1  —  exp  (~Yb/Yc) 
B  o  •  c 


[l-exp(-Yb/vc)] 


,PB,8PSK(vb^dYb  (5,22) 


where 


PB,8PSK'Yb'  =  §  Q[(6Yb)°-5sin(TT/8)] 
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Equal  Gain  Diversity 

This  diversity  technique  is  characterized  by  the  prop¬ 
erty  that  all  channels  have  exactly  the  same  gain  (3:1081). 
In  terms  of  Eq  (5.8) 

a^  =  1,  j  =  1,2,  ...  ,  L  (5.23) 

The  principle  of  this  method  is  that  all  the  received  chan¬ 
nel  signals,  fj(t),  are  added  together  and  the  resulting 
f(t)  is  formed  (See  Figure  5.3).  The  parameter  for  choosing 
the  optimum  signal  is  a  summation  of  the  signals  (i.e.,  sig¬ 
nal  strengths)  to  determine  which  symbol  signal  was  trans¬ 
mitted  . 


Figure  5.3  Equal  Gain  Diversity 


where 


a„  =  attenuation  factor  in  Kth  channel 

K 

=  energy  per  bit  in  the  transmitted  signal 
N„  =  noise  spectral  power  density  in  Kth  channel 
L  =  number  of  diversity  channels 
By  assuming  Rayleigh  fading  and  equal  noise  in  all  chan¬ 
nels,  the  p.d.f.  of  Yfc  for  equal-gain  diversity  combining  is 
developed  by  Schwartz  et  al.f  (20:456-7)  as 


where 


pD(V 


2L~1  LL  (Yb)L~1 

<2L'1>!  t  j<* 


(5.25) 


If  we  assume  all  diversity  channels  have  equal  mean  SNR 
over  the  short-term  fading,  ycK  =  Yc  then 


,L-1  T  L  .  \L-1 

,  .  2  L  (Yb}  (5.26) 

D  b  (2L-1) !  (YC)L 

Now  the  probability  of  bit  error  can  be  calculated  for 
the  equal-gain  diversity  combining  technique  using  DQPSK  and 
8-ary  PSK . 

For  equal-gain  diversity  using  DQPSK, 


Pn (DQPSK,  equal-gain  diversity)  =  f 
B  o 


2L_1  lL  (Yb)L_1 


(2L-1)!(YC)L 


‘PB,DQPSK^Yb^dYb  (5-27) 


where 


PB,DQPSK^Yb^  2Q[(4vbsin  (tt/8)  )  3 

For  equal-gain  diversity  using  8-ary  PSK, 

-  2L_1  LL  (y.  )] 


0.5- 


Pb(8PSK,  equal-gain  diversity)  =  f 


where 


o  (2L-1)!(yc)L 
’PB,8PSK^Yb)dYb 


PB, 8PSK ^Yb^  3  Q^6Yb)  sin(Tr/8)  3 


(5.28) 


Maximal  Ratio 

This  diversity  technique  is  similar  to  the  equal  gain 
method  as  the  received  channel  signals,  fj(t),  are  summed 
together  to  form  the  system  output  signal,  f(t);  however, 
unlike  the  equal  gain  technique,  maximal  ratio  requires  that 
the  gain  of  each  channel  is  weighted  by  a  factor  proportional 
to  its  received  signal  strength  (See  Figure  5.4).  Therefore, 
a  strong  received  signal  carries  a  larger  weight  and  con¬ 
tributes  more  to  the  system  output,  f(t),  than  a  weak  one. 
Each  received  channel  signal  is  multiplied  by  the  correspond¬ 
ing  complex -valued  (conjugate)  channel  gain  a^  exp  (j0K)  thus 
this  type  of  combiner  assumes  that  the  channel  attenuations 
and  phase  shifts  are  perfectly  known  (19:472). 

Maximal  ratio  diversity  is  one  of  the  most  commonly 
used  linear  diversity  combining  techniques  since  it  has  been 
shown  that  it  is  the  optimum  combining  method  with  respect 
to  probability  of  error  (3:1081  ,  9:2-8).  However,  to  obtain 
the  proper  weighting  factors  for  each  channel  requires  com¬ 
plex  instrumentation.  In  a  satellite  communications  down¬ 
link,  a  satellite  beacon  could  provide  the  information  needed 
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to  form  the  proper  weighting  factors.  For  this  reason,  some¬ 
times  equal  gain  or  selection  diversity  is  chosen  since  the 
additional  advantages  of  maximal  ratio  may  not  be  large 
enough  to  justify  the  additional  cost  of  the  equipment 
needed  to  achieve  the  necessary  maximal  ratio  combining. 


Analysis  of  Maximal  Diversity 

The  p.d.f.  of  the  output  SNR,  y^,  for  a  maximal  ratio 
combiner  over  short-term  fading  is  developed  by  Schwartz 
et  al . ,  (20:442-4).  By  taking  the  LaPlace  transform  of  the 
characteristic  function  of  yb  and  obtaining  the  inversion 
of  the  transform,  a  complicated  expression  for  PD(Yb)  is 
obtained  via  contour  integration.  However,  for  the  special 
case  when  all  average  channel  SNRs  are  equal  during  the 
short  term  fading,  this  expression  reduces  to 


(V 


L-l 


pD(Yb)  - - — f“  exp  <-Vb/vc)  (5.29) 

D  b  (L-l)!  (yc)L  b  c 

Now  the  probability  of  bit  error  can  be  calculated  for  the 
maximal-ratio  diversity  combining  technique  using  DQPSK  and 
8-ary  PSK . 

For  maximal  ratio  diversity  using  DQPSK, 

iL-1 

00  1  '■Yvy 

mavimal  v  a  4-  "i  /-»  —  P  -  - 

B 


oo  1  (Yb) 

PR(DQPSK,  maximal  ratio  diversity)  =  S  -  — - — —  exp(-yb/yc) 


o  (L-l)!  (yc)L 


■PB, DQPSK (yb)dyb  (5.30) 


where 


PB,DQPSK^Yb^  =  2Q^(4ybsin  [v/8)  ] 
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For  maximal  ratio  diversity  using  8-ary  PSK, 


Pb(8PSK,  maximal  ratio  diversity) 


f  \L-1 
j.°°_l _  (Yb) 

o  (L-1)!(yc)L 


■exp(-vb/Vc) 
•PB,8PSK(Yb)dYb  (5.31) 


where 

PB,8PSK(Yb)  =  |  Q[(6Yb)°,5sin(TT/8)] 


An  alternative  method  to  calculate  the  probability  of 
bit  error  for  maximal  ratio  diversity  systems  is  developed 
by  Proakis  in  Ref  19.  Proakis  (19:490)  develops  a  general 
result  for  the  probability  of  a  symbol  error  in  M-ary  PSK 
systems  using  maximal  ratio  diversity 


p  (-l)L-1d-u2)L 
K“  TT  (L-l)  5 


jL-1 
,  L-l 


b-u 


[  - 
L  M 


(M-l) 


u  sin(TT/M) _  -ucos  cos(tt/M)  s  -i  [  v 

r,  2  2,  a..  -i0 .5  cot  V,  2  2,  -i0 .5^  ’ 

Lb-u  cos  (tt/M)j  LL-u  cos  (rr/M)  J  ’ 

(5.32; 


where 

H  =  [yc/(l  +  Yc)D^*~*  for  coherent  PSK 
4  =  Yc/(1  +  Yc  for  DPSK 

E  \  aK2  f  £ 

Y  =  average  SNR  per  channel  =  - 


$  =  energy  in  transmitted  signal 


L  y 

Yb  =  average  bit  SNR  =  — 

K  =  log2  M 

The  bit  error  rate  for  DQPSK  with  maximal  ratio  diver¬ 
sity  is  derived  by  Proakis  assuming  a  Gray  code  (See  Figure 
5.5).  The  probability  of  bit  error  is  calculated  by 


P-n (DQPSK,  maximal  ratio  diversity)  =  Z  (K  errors)  »(P  )  (5. 

B  K=1  6K 


where 


33) 


=  probability  of  error  for  K  bit  errors 


The  resulting  expression  for  probability  of  bit  error 
for  DQPSK  vith  maximal  ratio  diversity  is 


Pg(DQPSK ,  maximal  ratio  diversity) 


0.5 


Cl- 


(2-p.2) 0  *5 


L-l 

£ 

K=0 


(5.34) 

4-2^i 

where 

H  =  Cyc/(1  +  Yc)]0,5 

Using  the  same  method  used  by  Proakis  for  obtaining 
Eg  (5.34),  the  probability  of  bit  error  for  an  8-ary  PSK 
signal  with  maximal  ratio  diversity  can  be  generated  by  eval¬ 
uating  Eg  (5.5)  where 


0 

/  2  P( 6) d0  is  defined  by  Proakis  [Eg  (7A.14),  Ref  19:532)] 


as 


S6 2  p(6)d0  =  ihl1  (_i — i  U.U-[(b/u2)-l1x2) 

J  2tt (L-l)  !  abL-1  b-^2  '  b0,5 


2\  ,-w2,0.5 


cot-1*  -  cot -x(— )  i  ); 

l-[b(^2)-l]x2)0,5  '  : 


b  =  1 


(5.35) 


where 


X,  = 


-u  cos  6 i 


(b-n2  cosei)0,5 


,  i  =  1,2 


0.5 


4  =  Cyc/(1  +  Yc)] 


Having  developed  the  equations  needed  to  calculate  the 
probability  of  bit  error  for  the  different  combinations  of 
diversity  techniques  and  digital  signaling  schemes,  it  was 


necessary  to  develop  a  computer  program  to  obtain  the  desired 
calculations.  In  order  to  reduce  run  time,  I  wrote  a  sepa¬ 
rate  program  in  Fortran  77  for  each  diversity  technique/dig¬ 
ital  signaling  scheme  combination.  Including  the  case  of  no 
diversity  (i.e.,  the  probability  of  bit  error  for  DQPSK  and 
8-ary  PSK  with  no  diversity) ,  there  are  seven  programs  used 
to  calculate  the  probability  of  bit  error  calculations. 

Choice  of  Range  of  Variables 

The  computer  results  were  to  be  plotted  as  the  probabil¬ 
ity  of  bit  error  versus  the  signal-to-noise  ratio  (SNR)  per 
bit.  Therefore,  a  range  for  the  SNR  per  bit  needed  to  be 
chosen.  Considering  the  range  (5-40  dB)  chosen  in  similar 
graphs  presented  in  Ref  19,  the  range  for  bit  SNR  was  chosen 
as  10-40  in  decibels.  A  determination  was  also  made  on  the 
maximum  number  of  diversity  channels  for  probability  of  error 
calculations.  Since  the  use  of  more  than  two  diversity  chan¬ 
nels  in  a  satellite  communications  diversity  system  is  rare, 

I  limited  the  probability  of  bit  error  calculations  for  the 
six  different  diversity /signaling  scheme  combinations  to  2, 

3,  and  4  diversity  channels.  Therefore,  the  final  results 
for  the  six  combinations  can  be  compared  over  the  bit  SNR 


range  of  10-40  decibels  for  2,  3,  and  4  diversity  terminals. 
Likewise,  these  results  can  be  compared  to  the  case  of  no 
diversity  to  determine  the  benefits  of  the  diversity  system. 


Equations  Used  to  Calculate  Probability  of  Bit  Error 

The  equations  needed  to  calculate  the  probability  of  bit 
error  for  all  of  the  diversity /signaling  scheme  combinations 
were  developed  in  Section  V.  The  equations  from  Section  V 
used  in  the  computer  programs  are  listed  below: 


DQPSK  w/no  diversity 
8-ary  PSK  w/no  diversity 
Selection  diversity  w/DQPSK 
Selection  diversity  w/8-ary  PSK 
Equal  Gain  diversity  w/DQPSK 
Equal  Gain  diversity  w/8-ary  PSK 
Maximal  Ratio  diversity  w/DQPSK 
Maximal  Ratio  diversity  w/8-ary  PSK 


Eq  (5.10) 
Eq  (5.11) 
Eq  (5.21) 
Eq  (5.22) 
Eq  (5.27) 
Eq  (5.28) 
Eq  (5.30) 
Eq  (5.31) 


All  of  the  equations  listed  above  require  the  solution 
to  an  indefinite  integral.  For  the  maximal  ratio  case,  an 
alternative  method  for  obtaining  the  probability  of  bit  error 
for  DQPSK  and  8-ary  PSK  is  stated  in  Eqs  (5.34)  and  (5.35), 
respectively.  The  difficulties  in  getting  the  required  par¬ 
tial  derivatives  for  this  method  resulted  in  the  deletion  of 


this  method. 


Solving  the  Indefinite  Integral 


The  choice  of  Fortran  77  as  the  programming  language 
for  the  programs  for  calculating  the  probability  of  bit  error 
was  primarily  due  to  the  mathematical  and  statistical  library 
available  in  Fortran.  In  order  to  solve  the  indefinite  inte¬ 
grals  needed  to  find  the  probabilities  of  bit  error,  I  used 
the  DCADRE  subroutine  from  the  IMSL  software  library.  The 
DCADRE  subroutine  performs  a  numerical  integration  of  a  func¬ 
tion  using  cautious  adaptive  Romberg  extrapolation. 

The  indefinite  integrals  in  the  equations  for  evaluating 

the  probability  of  error  were  approximated  by  using  a  piece- 

wise  summation  of  the  integration.  The  integral  from  zero 

to  infinity  is  approximated  by  the  piecewise  integral  from 
-15 

1  x  10  to  the  point  where  the  area  of  this  last  subinter¬ 
val  is  less  than  0.000001  of  the  summation  of  its  area  plus 
all  previously  calculated  subinterval  areas.  This  approxima¬ 
tion  works  satisfactorily  since  the  product  of  the  p.d.f.  of 
the  diversity  technique  and  the  conditional  probability  of 
bit  error  given  the  bit  SNR  is  always  positive  and  the  limit 
of  this  product  approaches  zero  as  bit  SNR  approaches  infinity 

The  DCADRE  subroutine  requires  that  the  integrand  of  the 
integral  be  supplied  by  the  user  in  the  form  of  a  function. 

For  each  diversity /signaling  scheme  combination,  the  product 
of  the  p.d.f.  and  conditional  probability  was  expressed  as  a 
Fortran  function  so  the  DCADRE  function  could  calculate  the 
subinterval  area.  In  order  to  calculate  the  conditional 


probability  of  a  bit  error  for  the  digital  signaling  scheme 
within  this  function,  another  subroutine  from  the  IMSL  soft¬ 
ware  library  was  needed. 


The  IMSL  subroutine  MERFD=DERF  was  used  to  evaluate  the 
mathematical  error  function  of  a  double  precision  argument. 
This  subroutine  computes  the  value  of  the  integral 

derf(Y)  =  (2/tt)  S*  exp(-t2)  dt  (6.1) 

o 


where 

Y  =  argument  of  the  mathematical  error  function 


By  performing  the  following  conversion,  the  Q  function 
defined  by  Eq  (2. 25. a)  can  be  expressed  as  a  function  of  the 
IMSL  subroutine  MERFD=DERF . 

Q (x)  =  1  -  [(derf (x/20,5)  +  l)/2]  (6.2) 

Therefore,  the  conditional  probabilities  in  the  integrand  of 
the  indefinite  integral  can  be  expressed  in  the  Q  function 
form  within  the  Fortran  function  needed  by  the  DCADRE  sub¬ 
routine  . 


Validation  of  Computer  Results 

The  computer  results  for  the  DQPSK  with  no  diversity 
and  DQPSK  with  maximal  ratio  diversity  cases  were  compared  to 
the  results  presented  by  Proakis  in  Ref  G,  Chapter  7.  For 
each  case,  the  results  were  similar,  although  my  probability 
of  bit  error  curves  are  slightly  shifted  to  the  right  of 
Proakis'  curves  (i.e.,  for  both  cases  my  results  show 


slightly  higher  probability  of  bit  error  rates  than  Proakis) . 
The  most  likely  cause  of  this  difference  is  the  approxima¬ 
tion  used  for  DQPSK  [Eq  (5.2)3.  However,  the  accuracy  of 
the  computer  results  is  well  within  the  range  to  validate 
the  approach  used  to  compute  the  probability  of  bit  error. 
Since  all  seven  programs  were  written  using  the  same  approach 
it  is  assumed  that  all  the  computer  results  are  reasonably 
accurate  and  in  particular,  the  results  are  accurate  for 
assessing  the  relative  performance  of  the  diversity  schemes. 


Computer  Program  for  Calculation  of  Rain  Attenuation 

In  addition  to  the  programs  for  calculating  the  proba¬ 
bility  of  bit  error,  a  program  was  developed  to  calculate  the 
attenuation  per  unit  distance.  A,  based  on  the  approximate 
relation  between  A  and  the  rain  rate  R,  where  A  =  aRb,  des¬ 
cribed  in  Eq  (3.2).  The  Values  for  the  coefficients  a  and 
b  were  obtained  for  three  different  rain  temperatures, 

-10°C,  0°C,  and  20°C.  These  values  were  obtained  from  Tables 
3.2,  3.3,  and  3.4  using  a  downlink  frequency  of  20  GHz  and  a 
Marshall-Palmer  (MP)  raindrop  size  distribution.  The  rain 
attenuation  per  unit  distance  will  be  calculated  as  decibels/ 
kilometers  and  the  rain  rate  will  vary  from  5  mm/hr  to 
100  mm/hr . 


VII .  Results 


The  data  obtained  from  the  computer  programs  described 
in  Section  VI  are  discussed  in  the  following  pages.  By  exam¬ 
ining  the  data,  an  evaluation  of  the  different  diversity/ 
signaling  scheme  combinations  can  be  obtained  to  determine 
the  benefits  of  using  a  diversity  network.  Comparisons  be¬ 
tween  the  different  diversity  techniques  will  illustrate 
which  technique  yields  the  maximum  benefits. 

Analysis  of  Diversity  Benefits 

By  using  the  HP-7229  plotter,  the  results  of  the  com¬ 
puter  programs  were  plotted  to  obtain  the  graphs  of  the  prob¬ 
ability  of  bit  error  versus  the  SNR  per  bit  for  each  diver¬ 
sity/signaling  scheme  combination.  These  graphs  are  depicted 
in  the  following  figures: 


DQPSK  and  8-ary  PSK 
w/No  Diversity 

Figure 

7.1 

Selection  Diversity 
w/DQPSK 

F igure 

7.2 

Selection  Diversity 
w/8-ary  PSK 

Figure 

7.3 

Equal  Gain  Diversity 
w/DQPSK 

Figure 

7.4 

Equal  Gain  Diversity 
w/8-ary  PSK 

Figure 

7.5 

Maximum  Ratio  Diversity 
w/DQPSK 

Figure 

7.6 

Maximum  Ratio  Diversity 
w/8-ary  PSK 

Figure 

7.7 
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By  looking  at  Figures  7.1  -  7.7,  it  is  obvious  that 
the  use  of  a  diversity  technique  greatly  reduces  the  proba¬ 
bility  of  bit  error  in  a  Rayleigh  fading  channel.  It  is 
also  apparent  that  regardless  of  which  diversity /signaling 
scheme  is  used,  the  probability  of  bit  error  reduces  as  the 
order  of  diversity  of  channels  increases.  A  more  indepth 
look  at  the  diversity  benefits  are  depicted  in  Tables  7.1  - 
7.4. 

As  displayed  in  Tables  7.1  and  7.2,  the  benefits  of  di¬ 
versity  techniques  are  significant.  Table  7.1  shows  the 
benefits  of  the  selected  diversity  techniques  with  DQPSK  com¬ 
pared  to  DQPSK  without  diversity  and  Table  7.2  shows  the 
benefits  of  the  selected  diversity  techniques  with  8-ary  PSK 
compared  to  8-ary  PSK  without  diversity.  For  each  case,  as 
the  average  SNR  per  bit  increases,  the  benefits  of  the  diver¬ 
sity  network  increases.  For  example,  given  a  received  aver¬ 
age  bit  SNR  of  40  dB  for  a  DQPSK  (8-ary  PSK)  w/no  diversity 
system,  the  minimum  diversity  gain  is  18.25  dB  (20.0  dB) 
when  using  selection  diversity  with  two  channels  and  the  max¬ 
imum  diversity  gain  is  26  dB  (27.25  dB)  when  using  maximal 
ratio  diversity  with  four  channels.  Likewise,  for  a  received 
average  bit  SNR  of  60  dB,  the  diversity  gain  for  selection 
diversity  with  two  channels  is  28.25  dB  (30.0  dB)  and  the 
diversity  gain  for  maximal  ratio  diversity  with  four  channels 
is  40.5  dB  (41 .75  dB) . 

The  data  displayed  in  Tables  7.3  and  7.4  allows  for  a 
comparison  of  all  six  diversity/signaling  scheme  combinations 


Two  significant  conclusions  can  be  made  from  this  datai 


1.  The  maximal  ratio  diversity  technique  provides  the 
most  diversity  gain,  followed  by  equal  gain  diversity  and 
lastly  selection  diversity. 

2.  The  diversity/signaling  scheme  combinations  with 
8-ary  PSK  outperform  the  same  combinations  with  DQPSK . 

(The  latter  conclusion  may  partially  result  from  the  approx¬ 
imation  used  for  DQPSK) .  Based  on  these  results,  the  optimum 
diversity /signaling  scheme  combination  of  the  six  is  the  max¬ 
imal  ratio  diversity  technique  with  8-ary  PSK. 

Other  important  considerations  regarding  the  selection 
of  a  diversity  technique  for  a  satellite  communications  sys¬ 
tem  is  the  cost  of  implementing  the  diversity  technique.  As 
denoted  in  Section  V,  maximal  ratio  diversity  requires  the 
use  of  a  satellite  beacon  to  properly  weigh  the  received  chan 
nel  signals.  Although  maximal  ratio  diversity  provides  the 
best  diversity  gain,  it  is  also  the  most  costly  to  imple¬ 
ment  due  to  the  additional  equipment  needed  to  obtain  the 
satellite  beacon  information  and  translate  it  into  correct 
weighing  factors.  Thus,  there  is  a  direct  trade-off  between 
cost  and  performance  in  the  selection  of  the  diversity  tech¬ 
nique  and  the  number  of  earth  stations  used  in  the  satellite 
communications  diversity  system. 

Analysis  of  Rain  Attenuation  for  Downlink 

The  rain  attenuation  for  a  satellite  communications  down 
link  can  be  estimated  from  the  data  presented  in  Figures  7.8 


and  7.9.  In  Figure  7.8,  the  rain  attenuation  per  unit  dis¬ 
tance  (dB/km)  can  be  approximated  by  knowing  the  rain  rate 
for  rain  temperatures  of  20 °C,  0°C,  and  -10°C.  Given  the 
rain  rate  and  the  angle  of  elevation  of  the  receiving  satel¬ 
lite  terminal,  the  effective  path  length  (in  km)  can  be  esti¬ 
mated  using  Figure  7.9. 

Therefore,  the  approximate  rain  attenuation  (in  dB)  can 
be  obtained  by  multiplying  the  effective  path  length  deter¬ 
mined  from  Figure  7.9  and  the  rain  attenuation  per  unit  dis¬ 
tance  indicated  from  Figure  7.8.  By  knowing  the  approximate 
rain  attenuation,  its  effect  of  the  probability  of  a  bit 
error  can  be  determined  from  Figures  7.1  -  7.7.  By  subtrac¬ 
ting  the  rain  attenuation  from  a  given  received  bit  SNR  for 
a  diversity /signaling  scheme  for  a  satellite  communications 
downlink  without  rainfall  effects,  the  corresponding  proba¬ 
bility  of  bit  error  can  be  obtained.  The  increase  in  the 
probability  of  bit  error  is  the  result  of  the  rain  attenua¬ 
tion.  Therefore,  the  effects  of  rain  attenuation  for  differ 
ent  rain  rates  and  different  elevation  angles  on  the  six  dif 
ferent  diversity /signaling  scheme  combinations  can  be  esti¬ 
mated  using  Figures  7.1  -  7.9. 

For  example,  assume  a  maximal  ratio  diversity  system 
with  two  channels  using  8-ary  PSK  for  a  downlink  earth  sta¬ 
tion  with  an  angle  of  elevation  of  40° operating  during  clear 
weather  conditions.  For  an  average  received  SNR  per  bit  of 
20  dB,  the  probability  of  bit  error  obtained  from  Figure  7.7 
is  2.807  x  10“5.  Now  assume  a  rain  rate  of  25  mm/hr  and  a 


rain  temperature  of  20°C.  From  Figure  7.8,  the  rain  attenua 
tion  is  ~2.6  dB/km  and  from  the  effective  path  length  of 
~4.4  km  is  obtained  from  Figure  7.9.  Thus,  the  rain  atten¬ 
uation  is  ~11.44  dB  and  the  effective  average  received  SNR 

per  bit  for  the  diversity  system  is  now  ~9.5  dB  with  a  prob- 

-3 

ability  of  bit  error  of  2.786  x  10  .  Obviously,  the  detri¬ 

mental  effect  of  rain  attenuation  on  the  performance  of  sat¬ 
ellite  communications  systems  is  very  significant. 


Figure  7.1  DQPSK  and  8-ary  PSK  w/No  Diversity 
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Average  tt=  Needed  to  Achieve  Prob.  of  Bit  Error  Performance 


Comparison  of  Selected  Diversity  Techniques  w/8-ary  PSK  to  8-ary  PSK  w/fao  Diversity 


Comparison  of  Selected  Diversity  Techniques  v/DQPSK 
with  Respect  to  Probability  of  Bit  Error 
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VIII .  Conclusions  and  Recommendations 

Conclusions 

On  the  basis  of  the  computer  results  of  the  performance 
of  diversity  techniques  in  the  downlink  reception  equipment 
for  satellite  communications  systems,  the  following  conclu¬ 
sions  are  made: 

1.  The  best  link  parameter  to  use  in  the  process  of 
determining  the  received  signal  in  a  diversity  system  is  the 
received  channel  SNR  or  some  variation  of  it  such  as  average 
SNR  per  bit. 

2.  Of  the  three  diversity  techniques  evaluated,  maximal 
ratio  diversity  provided  the  optimum  diversity  benefits  fol¬ 
lowed  by  equal  gain  diversity  and  lastly  selection  diversity. 
However,  all  three  diversity  methods  provide  significant 
improvement  of  bit  error  performance  when  compared  to  a  sim¬ 
ilar  system  without  diversity.  These  diversity  combining 
techniques  can  be  used  to  combat  rain  attenuation  or  to 
improve  the  overall  reception  characteristics. 

3.  Regardless  of  which  diversity  technique  is  used,  as 
the  number  of  channels  used  increases,  tne  bit  error  perform¬ 
ance  improves . 

4.  Maximal  ratio  diversity  requires  the  use  of  a  satel¬ 
lite  beacon  while  equal  gain  and  selection  diversity  do  not. 
Thus,  additional  equipment  is  needed  in  a  maximal  ratio 
diversity  reception  system  thereby  increasing  the  cost. 


5.  The  diversity /signaling  combinations  using  8-ary  PSK 
provided  a  slightly  better  bit  error  performance  than  its 
counterpart  using  DQPSK . 

6.  The  effects  of  rain  attenuation  on  a  diversity 
downlink  reception  system  can  be  estimated  by  applying  the 
results  of  this  study. 

Recommendations 

Based  on  the  observations  made  during  this  investiga¬ 
tion,  the  following  recommendations  are  proposed  for  further 
study: 

1.  This  study  assumed  perfect  synchronization;  the  ef¬ 
fects  of  the  loss  of  synchronization  on  bit  error  performance 
should  be  investigated. 

2.  The  results  of  this  study  could  be  applied  to  rain¬ 
fall  statistics  to  obtain  theoretical  long-range  uptime  and 
downtime  statistics  for  downlink  diversity  reception  systems 
located  in  various  geographical  locations. 
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